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ABSTRACT 


An  experimental  study  with  two  specific  aims  was  conducted  at  the  Air  Force 
Weapons  Laboratory:  (1)  to  show  the  relationship  between  temperature  and 
velocity  of  longitudinal  ultrasonic  waves  through  several  plastics  and  noe  - 
cone  materials  and  (2)  to  establish  a  convenient  apparatus  fcr  the  determina¬ 
tion  of  elastic  properties  of  msterialB.  The  velocity-temperature  data  were 
needed  for  use  in  another  research  project  where  certain  Inputs  were  necessary 
to  more  completely  describe  materials  under  shock-loaded  conditions.  The 
experimental  procedure  was  based  on  measuring  the  transit  times  requir  d  for 
ultrasonic  waves  in  the  low-megacycle  range  to  pass  through  samples  01  differ¬ 
ent  thicknesses.  The  resulting  accuracy  for  most  of  the  measurements  is 
within  about  1  to  2  percent.  The  temperature  range  of  the  measurements  extend¬ 
ed  between  room  temperature  and  about  125aC,  with  a  few  measurement?  to  250°C. 
Velocity-temperature  curves  were  determined  for  nylon,  low-  and  high-density 
polyethylene,  plexiglass  (polymethylmethacrylate),  Delrin  Acetal,  and  teflon. 
Curves  were  also  drawn  for  several  re-entry  vehicle  nose-cone  materials,  in¬ 
cluding  chopped  nylon  phenolic.  Castable  124,  Avcoat  1  and  19,  pyrolytic 
graphite.  General  Electric  Phenolic  Fibre  Jlass,  t'henolic  Carbon,  Tape  Wound 
Nylon  Phenolic,  Rad  S8B,  Rad  60,  and  some  solid  epoxy  foams.  Less  complete 
data  were  obtained  for  Avco  phenolic  fiber  glass  and  Oblique  Tape  Wound  Refrasil. 
The  angular  dependence  of  velocity  was  measured  for  several  layered  materials, 
and  it  was  found  that  the  velocity  was  very  dependent  on  the  direction  of  propa¬ 
gation.  For  most  of  the  materials  studied  the  velocity  was  found  to  decrease 
more  or  less  linearly  with  Increasing  temperature  and  normally  was  abouc  9  - 
20  percent  low sr  then  at  room  temperature. 
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ABBREVIATIONS  AND  SYMBOLS 

Young's  modulus,  ratio  of  applied  stress  to  fractional  extension. 
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Elastic  stiffness  constants  which  are  the  coefficients  of  the  strains  in 

the  generalized  form  of  Hooke's  law.  In  matrix  form  o  *  [c  u,,. 

ij  I  ill  j-J 

Components  of  strain.  \  / 

Components  of  stress. 


Lafie  constant  for  isotropic  materials.  Identical  with  the  shear  modulus  G. 
Laifte  constant  appearing  in  the  stres6~strain  matrix  for  isotropic  materials. 


Longitudinal  velocity,  defined  in  isotropic  media  as  c^ 

Shear  velocity,  defined  in  isotropic  media  as  cfc 

/3Po\  " 

Hydrodynamic  sound  velocity,  defined  as  c£  ■l-jjjj-'l  ,  wh 
and  Pq  is  the  pressure  at  normal  conditions. 


X  +  2V 


where  s  is  entropy 


Shock  velocity  defined  as  Ug  •!— I,  along  the  actual  path  in  the  pressure- 
density  plane.  VD/ 

Acoustic  impedance,  defined  as  product  of  velocity  (longitudinal  or  shear) 
and  initial  sample  density  p. 

Particle  velocity. 

Polason'a  ratio,  ratio  of  latrrnl  contraction  to  longitudinal  extension 

of  a  specimen . 

Specific  heat  at  constant  pressure. 

Attenuation  of  an  acoustic  wave  in  uepers.  /  v 
Volume  coefficient  of  thermal  expansion  B  -  yrfjL* 


Cruenelaen  ratio,  defined  as  y  -  *#her*  P  **  the  PreMure*  E  th* 

internal  energy,  and  V  is  specific  volume. 

o  -  o 


Compression  defined  as  u '  - 


-,  where  is  the  initial  sample  density 


and  o  is  the  density  behind  the  pressure  wave. 
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SECTION  I 
INTRODUCTION 

The  first  fact  which  becomes  apparent  to  a  researcher  desiring  information 
on  the  velocity  of  ultrasound  penetrating  through  a  particular  material  is  that 
there  is  a  shortage  of  publish  1  specific  parametric  data  on  the  phenomenon. 

Such  parameters  as  the  frequency  of  the  acoustic  waves  used  and  the  temperature 
at  which  measurements  were  made  are  noticeably  lacking.  For  plastics  and  re¬ 
lated  materials,  in  particular,  there  is  a  scarcity  of  information  indicating 
the  relationship  between  velocities  and  temperatures. 

For  this  reason  an  investigation  was  conducted  at  the  Air  Force  Weapons 
Laboratory  to  determine  as  functions  of  temperature,  the  velocities  of  longi¬ 
tudinal  ultrasonic  waves  through  several  plastic  materials,  including  some  of 
the  nonmetallic  materials  used  in  the  nose  cones  of  ballistic  reentry  vehicles. 
This  information  was  needed  as  input  data  in  other  research  programs  which  were 
concerned  with  the  dynamic  response  of  materials  t  5  shock  loading.  Without 
published  data  showing  the  effect  of  temperature  upon  velocity  one  might  assume 
the  velocitv  to  be  constant  over  fairly  large  temperature  ranges.  This  is 
indeed  a  very  broad  assumption  which  the  results  of  this  study  have  shown  to  be 
invalid. 

In  view  of  the  accuracy  of  other  input  data  and  the  reproducibility  and 
inhomcgt.-.eity  of  sample  material  it  was  determined  that  a  precision  of  0.5Z  iv. 
velocity  .slues  was  needed.  Considering  this  and  the  availability  of  equipment 
it  was  decided  to  use  a  simple  and  straightforward  direct  timing  method  in  which 
transit  ?iraes  through  samples  of  different  thicknesses  are  measured. 

Before  proceeding  further,  a  word  on  definitions  is  in  order.  The  term 
"ultrasonic"  is  »'ow  commonly  used  to  describe  vibrational  waves  with  frequencies 
higher  than  th'-'se  normally  audible  to  the  ear,  i.e.,  any  acoustic  wave  of  a 
frequency  higher  than  approximately  20  kc  is  said  to  be  ultrasonic.  In  some 
older  works  the  word  "supersonic"  was  used  to  mean  the  same  thing,  but  such 
usage  is  no  longer  sanctioned  by  convention.  Similarly,  "ultrasonics"  is  de¬ 
fined  as  the  te.  .mology  of  sound  at  frequencies  above  the  audio  range,  and 
"supersonlcs"  is  now  defined  to  be  the  general  subject  covering  phenomena 
associated  with  speed  higher  than  the  speed  of  sound  (ss  in  the  case  of  aircraft, 
projectiles,  etc).  This  study  is  concerned  only  with  "ultrasonics." 
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SECTION  II 
EXPERIMENTAL  METHOD 

The  technique  of  determing  velocity  by  measuring  the  time  required  for 
pulses  to  traverse  samples  01  differ*,  t  thicknesses  was  chosen,  whenever  possi¬ 
ble,  to  eliminate  complicating  considerations  such  as  transducer  and  circuitry 
delay  times,  and  interfacial  effects. 

The  resulting  sample  arrangement  is  best  explained  by  referring  to  figure 
1.  If  the  transit  time  form  A  to  B  through  the  thick  sample  is  tj,  and  the 
transit  time  through  the  thin  sample  is  t2,  the  velocity  of  the  pulse  through 
the  sample  material  can  be  calculated  as, 

c  .  “  .  h.  a) 

ci  At  tj  -  t2 

where  c£  is  the  longitudinal  bulk  velocity.  If  the  same  aluminum  block*  and 
the  same  transducers  arc  used  while  Baking  the  measurements  of  tj  and  t2,  it 
is  not  necessary  to  determine  the  transducer  delay  times  or  the  transit  time 
through  the  aluminum  buffers,  since  these  constant  delays  disappear  on  comput¬ 
ing  the  difference  in  times,  tj  -  t2.  This  method,  of  course,  assumes  that  the 
acoustic  units  can  be  reproducibiy  assembled  and  that  the  difference  in  transit 
times  is  due  only  to  the  difference  in  sample  thicknesses. 

Figure  2  shows  an  Improved  bonding  clamp  based  on  Sullivan's  design  (Ref.  1) 
which  was  devise*  to  achieve  this  reproducibility.  It  was  found  that  with  this 
technique  the  deviation  in  time  due  to  assembly  was  generally  within  +  O.Olu  ~ 
sec.  The  transmission  blocks  were  cut  from  2  inch  diameter  bar  stock  of  2024 
aluminum  alloy.  The  ends  of  the  rods  were  machined  and  polished  so  as  to  be 
flat,  parallel,  and  smooth,  thereby  reducing  diffraction  and  edge  effects.  The 
length  of  the  rods  used  in  this  work  were  varied  to  suit  the  individual  experi- 
ment  and  were  normally  on  the  order  of  3/4  to  2  Inches  long. 

*The  function  of  the  alusdnum  block  is  to  make  the  overall  transit  time  greater 
than  the  Juration  of  the  ultrasonic  pulse  so  that  electrical  noise  generated 
by  ':he  transmitted  pulse  does  not  interfere  with  the  received  pulse. 


2 


Figure  1.  Diagrams  of  Typical  Acoustic  Assembly 
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The  samples  were  also  cut  as  right  circular  cylinders  of  approximately 
2-inch  diameter.  The  exact  shape  of  the  sample  is  immaterial,  provided  that  the 
lateral  and  longitudinal  dimensions  are  large  enough  that  the  assumption  of  an 
infinite  medium  is  valid.  Tu,  Brennan,  and  Sauer  (Ref.  2)  have  found  ultrasonic 
velocities  to  be  independent  of  sample  size  provided  the  radius  of  the  sample  is 
at  least  2.5  times  as  great  as  the  wavelength  of  the  stress  disturbance.  The 
samples  and  transmission  blocks  used  in  the  present  study  have  radii  that  are 
generally  six  or  more  times  greater  than  the  wavelengths  used.  This  ensures 
the  one  dimensionality  of  the  hydrodynamic  wave  so  that  the  longitudinal  bulk 
velocity  is  essentially  the  quantity  measured. 

For  the  required  accuracy  it  was  imperative  that  the  two  end  faces  of  the 
sample  be  flat  and  parallel  to  reduce  diffraction  and  dispersion  effects. 

This  was  achieved  with  varying  degrees  of  success  depending  on  the  particular 
sample  material  and  the  technique  used  in  preparing  it.  Normally,  however,  it 
was  found  that  the  variation  of  sample  thickness  was  so  small  that  the  contri¬ 
bution  of  the  resulting  error  to  the  overall  error  could  be  neglected.  This 
source  of  error,  which  contributed  at  most  a  deviation  of  0.5  percent,  will  be 
treated  in  more  detail  in  a  later  section. 

The  transducers  used  to  produce  the  acoustic  vibrations  are  composed  of 
lead  zirconate/titanate  mixtures  f?-med  in  the  shape  of  circular  disks  and 
designed  to  operate  in  the  thickness  mode.  Their  thicknesses  varied  depending 
upon  the  desired  resonant  frequency,  e.g. ,  ^he  thickness  corresponding  to  a 
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C i 026  iaCuca *  In  most  cases, 
the  transducers  were  ope  ated  at  the  fundamental  frequencies  of  1.5,  3.0,  5.0 
and  10.0  megacycles  per  second,  or  the  respective  harmonics. 

For  most  of  the  measurements  a  light  transformer  oil*  was  used  as  the 
coupling  agent  between  the  various  components  of  the  acoustic  units.  For 
temperature  measurements  (to  temperatures  of  about  125*0  the  same  oil  was  used 
as  the  liquid  la  a  thermostatically  controlled  bath  used  to  control  sample  tem¬ 
peratures.  For  temperatures  higher  than  this  and  for  samples  which  were  oil 
absorbent,  a  thermostatically  controlled  oven  was  used  for  the  temperature 
measurements. 


*"Special  Marco 1  42-46",  produced  by  Humble  Oil  and  Refining  Company. 
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The  acoustic  unit  shown  in  figure  2  is  composed  of  two  ceramic  transducers 
(only  the  upper  one  is  visible),  a  sample,  and  an  aluminum  transmission  block. 
The  cylindrical  aluminum  electrodes  of  the  clamp  are  mounted  on  ball  joints 
which  ensures  that  the  components  of  the  acoustic  unit  are  held  together, 
accurately  aligned,  and  under  constant  spring  pressure. 

For  many  of  the  experiments  only  one  buffer  rod  was  used  in  calculating 
velocities.  A  grounding  electrode  was  therefore  provided  for  the  side  of  the 
transducer  next  to  the  sample,  since  the  plastic  samples  are  nonconducting. 

This  was  a  thin  (0.0017  inch  thick)  piece  of  aluminum  foil  which  can  be  seen  in 
figure  2  between  the  sample  and  the  upper  transducer.  The  wire  contact  was 
pressed  firmly  against  the  foil  with  spring  pressure  to  maintain  an  adequate 
ground  lead.  For  samples  which  were  highly  attenuating  or  for  which  it  was  not 
possible  to  obtain  a  sufficiently  thick  sample,  the  foil  was  eliminated  and  two 
aluminum  buffer  rods  were  used  instead  of  one.  This  modification  required  the 
transit  time  measurement  of  only  one  sample,  and  the  velocity  was  computed  as 
in  equation  (1)  with  Xj  equal  to  zero,  tj,  the  transit  time  with  the  sample 
Inserted,  and  t2,  the  transit  time  through  the  assembly  without  the  sample. 
Comparison  between  the  two  methods  was  good,  providing  proper  precaution  was 
taken  in  assembling  the  various  acoustic  units. 

The  main  difference  in  velocity  measurements  which  could  result  from  the 
use  of  these  tv  o  different  methods  is  that  occurring  because  of  the  difference 
in  the  number  of  sample-buffer  interfaces  with  and  without  the  sample  inserted 
for  the  latter  method.  However,  this  effect  appeared  to  be  insignificant  in 
comparison  to  the  Inherent  limitations  of  the  technique,  and  cross  checks 
between  the  two  methods  on  Identical  samples  agreed  regularly  to  wlthing  IX. 

To  distinguish  between  the  various  methods  referrsd  to  in  the  text,  the  tech¬ 
nique  of  using  different  samples  to  calculate  the  velocity  will  be  referred  to 
as  the  "comparison  method",  and  that  of  using  only  one  sample  with  the  corres¬ 
ponding  system  delay  times  will  be  referred  to  as  the  "single  sample  technique". 

Since  the  major  objective  of  the  experiment  was  to  determine  the  longi¬ 
tudinal  velocity  as  a  function  of  temperature,  it  was  necessary  to  know  the 
heating  rate  of  the  samples.  This  was  a  very  important  consideration,  since 
in  order  to  measure  the  true  velocity-temperature  dependence  it  was  necessary 
that  all  temperature  gradients  vanished  along  the  path  of  the  acoustic  beam. 
With  this  in  mind,  a  study  of  the  temperature-time  dependence  of  the  samples 
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was  conducted  to  determine  the  time  required  to  reach  temperature  equilibrium 
in  the  sample.  The  detailed  results  are  presented  in  Appendix  I,  but  a  few 
words  about  the  temperature  measurements  will  be  given  here. 

Sample  temperatures  were  determined  by  reading  the  bath  or  oven  tempera¬ 
tures  to  the  nearest  degree  by  means  of  a  mercury  thermometer  or  chromel-alumel 
thermocouple,  respectively.  The  agreement  between  the  two  methods  of  tempera¬ 
ture  measurement  was  found  to  be  within  18C  for  the  reported  temperature  range. 
The  results  of  the  work  shown  in  Appendix  I  indicated  that  approximately  one 
hour  was  necessary  before  the  sample  reached  thermal  equilibrium.  The  thinner 
samples  were  allowed  somewhat  less  time  for  the  temperature  to  stabilize,  but 
in  no  case  was  less  than  one  half  hour  allotted  for  temperature  stabilization. 
With  these  precautions  taken,  it  is  estimated  that  the  temperature  of  the 
sample  was  uniform  to  within  1°C  and  to  within  1*C  of  the  desired  environment. 

The  ultrasonic  pulse  generating  equipment  is  an  Aienberg*  pulsed  oscillator 
wiiich  generates  pulses  a  few  microseconds  long  with  a  carrier  frequency  in  the 
low  megacycle  range.  The  attenuators  indicated  in  figure  3  are  also  made  by 
Arenberg  and  are  designed  to  give  from  1  to  122  db  of  attenuation  in  1  db  in¬ 
crements.  Their  main  use  is  to  prevent  overloading  of  the  oscilloscope  for 
constant  input  signals  and  to  allow  relative  attenuation  measurements  to  be 
made.  The  time-mark  generator**  produces  1  microsecond  markers  and  is  cali¬ 
brated  to  provide  an  error  of  less  than  one  part  per  million.  The  oscilloscope 
used  for  all  measurements  is  a  well-adjusted  dual  beam  oscilloscope  with  two 
type  CA  dual  trace  plug-in  units.#  The  Arenberg  pulse  oscillator  was  rared  for 
an  output  of  600  volts,  peak-to-peak,  so  that  most  measurements  were  possible 
without  the  use  of  additional  amplifiers.  When  additional  amplification  was 
necessary,  the  amplifier  was  Inserted  on  the  receiving  branch  between  the 
attenuator  and  the  receiving  transducer. 

In  many  of  the  pulse  methods  of  measuring  ultrasonic  velocities  which  have 
been  reported  in  the  literature,  the  received  pulses  are  rectified  and  the  pulse 


*Arenberg  Ultrasonic  Laboratory,  Model  PG-650C,  Boston,  Massachusetts. 

**Hickok  Model  1817. 
iTektronic  Type  555. 
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envelope  observed.  Eros  and  Kelts  (Kef.  3)  have  described  a  single-transducer 
method  baaed  on  the  observation  of  the  unrectlfled  pulses.  This  vas  the  practice 
used  In  the  present  work.  In  making  velocity  determinations  by  such  a  timing 
method  it  Is  necessary  to  have,  on  both  the  transmitted  and  received  pulses, 
reference  points  between  which  the  elapsed  time  can  be  measured.  Therefore,  the 
unrectlfled  pulses  were  observed  in  their  entirety  on  the  oscilloscope,  so  that 
the  time  could  be  measured  from  one  particular  oscillation  of  the  transmitted 
pulse  to  the  corresponding  oscillation  of  the  received  pulse.  During  the  course 
of  this  investigation,  several  different  time-measuring  techniques  were  used, 
and  a  method  was  developed  which  geve  the  best  reproducibility  and  reliability 
and  for  which  most  of  the  data  are  reported. 

In  the  resulting  technique,  ultrasonic  pulses  approximately  3  u  sec  in 
duration  with  a  carrier  frequency  in  the  lew  megacycle  range  were  used  to  measure 
velocities.  The  received  and  transmitted  signals  were  delayed  on  the  oscillo¬ 
scope  by  a  calibrated  delay  sweep  potentiometer  in  order  to  identify  two 
corresponding  peaksvnear  the  middle  of  the  respective  pulses  as  references. 

The  transit  time  was  t'r&n  determined  by  counting  and  extrapolating  the  number 
of  time  marks  between  the  two  reference  peaks.  Figure  4  illuatrates  the  method. 

The  illustration  shows  three  hypothetical  oscilloscope  traces  with  a  fixed 
time  relationship  to  each  other.  They  represent,  from  top  to  bottom,  the  trans¬ 
mitted  pulse,  the  received  pulse,  and  a  series  of  1  u  sec  time  marks.  The 
desired  time  interval,  t,  la  that  between  the  4th  positive  peaks  of  the  trans¬ 
mitted  and  received  pulses.  This  tins  is  obtained  by  taking  readings  of  the 
potentiometer  corresponding  to  the  six  events  marked  p1  through  p6  *r>d  calcu¬ 
lating  the  elapsed  time  by  assuming  that  tbs  oscilloscope  delay  circuit  is 
sufficiently  linear  over  any  1  u  sec  interval  to  permit  linear  interpolation. 

With  this  time  measuring  techniqus,  successive  measurements ,  using  the  asms 
acoustic  unit,  regularly  agreed  to  within  +  0.004  ii  sec,  mad  most  measurements 
agreed  even  closer.  Host  of  the  velocities  reported  here  were  determined  by 
this  technique. 

In  certain  eases  the  thicker  samples  of  some  of  the  materials  studied  caused 
either  high  attenuation  or  excessive  distortion  of  the  acoustic  waves.  When 
this  condition  was  encountered,  the  velocity  measurement  method  empl  ed  by 
lolls  and  Slack  (tef.  4),  or  the  single  sample  method  described  earlier  was  used 
to  velocity  determinations  using  only  one  thin  sample  of  the  material.  In 
the  echo  method  of  Voile,  a  sample  was  placed  between  two  transmission  blocks 
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Figure  4.  Schematic  Representation  Illustrating  an  Interpolation 

Tine-Measuring  Technique 
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with  a  transducer  at  each  end  of  the  assembly,  and  the  transit  time  for  a  pulse 
to  pass  from  end  to  end  of  the  assembly  was  measured.  Then  the  time  required 
for  an  echo  to  return  from  the  sample  interface  to  the  sending  transducer  was 
measured  for  each  of  the  blocks  by  transmitting  from  each  of  the  transducers  in 
turn.  The  average  of  the  block  transmission  times  was  then  subtracted  from  the 
end-to-end  transit  time  to  obtain  the  transit  time  through  the  sample.  This 
method  was  infrequently  used  and  the  results  obtained  by  it  were  used  only  for 
comparison  purposes. 
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SECTION  III 
EXPERIMENTAL  RESULTS 

1*  Longitudinal  Velocity  In  Ccmaaon  Plastics 

The  first  common  plastic  studied  over  a  range  of  temperatures  was  a  nylon.* 
Temperatures  ranged  from  a  low  of  25°C  to  a  high  of  114eC  and  the  resultant 
velocity  values  are  shown  in  figure  5.  The  frequency  used  fo^  most  of  the 
measurements  was  3.0  Me,  although  velocity  checks  at  frequencies  of  1.5,  4.33 
8nd  10.0  Me  agreed  within  experimental  variation  with  the  values  at  3.0  Me,  thus 
Indicating  that  the  material  is  not  dispersive  in  this  frequency  range. 

Another  nylon  (Delrin  Acetal)  was  also  analyzed  as  a  function  of  temperature, 
and  the  results  are  plotted  in  figure  6.  The  density  for  this  material  uas 
1.43  g/cc,  and  no  apparent  dependence  with  frequency  observed  for  the 
longitudinal  velocity.  In  anticipation  of  velocity-density  measurements  re¬ 
ported  in  a  later  section,  it  is  of  interest  to  point  out  the  differences  in 
velocity  for  the  two  nylon  materials.  The  higher  density  sample  exhibits  a 
lower  value  of  longitudinal  velocity  at  room  temperature  than  the  lower  density 
material.  For  materials  which  are  Identical  in  composition  and  vary  only  in 
density  characteristics,  the  effect  is  generally  the  opposite,  which  is  shown 
by  the  following  two  materials  and  by  the  results  to  be  presented  later. 

Two  different  density  polyethylene  samples  were  analyzed  and  the  results 
are  shown  in  figures  7  and  8.  It  is  noted  that  the  velocity  at  room  temperature 
for  the  0.964  g/cc  (Fig.  8)  sample  Is  greater  than  that  for  the  0.916  g/cc 
sample,  agreeing  with  measurements  obtained  In  velocity-density  studies  of 
polyethylene  (Kef.  5).  The  temperature  range  used  here  was  necessarily  limited 
to  about  70*C  to  avoid  smiting  or  distorting  the  samples. 

The  longitudinal  wave  velocity  in  Plexiglas  (polywethylr Mthacrylate,  density 
1.19  g/cc)  <«  plottad  in  figure  9.  As  indicated  in  the  previous  graphs,  simi¬ 
lar  velocities  are  obtained  et  various  frequencies  studied  so  that,  again, 
dispersion  is  not  evident,  within  the  experimental  errors  present. 

The  lest  common  plastic  investigated  was  Teflon  (density  2.19  g/cc).  For 
this  material  the  temperature  range  was  extended  to  below  room  temperature  in 


*Poly-Panco  Nylon  (density  1.15  g/cc) 
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Figure  7.  Longitudinal  Wave  Velocity  ve  Temperature  in  Low  Density 
Polyethylene  (density  0.916  g/cc) 
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Longitudinal  Wave  Velocity  vs  Temperature  in  High  Density  Polyethylene 

(density  0.964  g/cc) 


Figure  9.  Longitudinal  Wave  Velocity  vs  Temperature  in  Plexiglass 
(polymethylmethacrylate,  density  1.19  g/cc) 
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order  to  obtain  a  more  accurate  picture  of  the  shape  of  the  velocity  curve  in 
the  range  between  20  and  30°C.  As  can  be  seen  in  figure  10,  there  is  a  fairly 
sharp  inflection  in  the  velocity  curve  in  this  region.  It  is  "oted  that  Teflon 
exhibits  a  well  known  second  order  transition  in  the  molecular  structure  at  a 
temperature  of  20°C  (Ref.  6),  so  that  it  is  to  be  expected  that  a  change  in  the 
mechanical  properties  (and  consequently  the  velocity)  would  be  observed  as  the 
temperature  is  varied  over  the  transition  region.  The  velocities  measured  in 
Teflon  appear  to  be  somewhat  more  frequency  sensitive  than  for  the  previously 
listed  materials.  However,  the  difference  is  not  sufficient  to  justify  separ¬ 
ate  plots ,  and  all  measurements  for  the  three  frequencies  used  are  plotted  on 
the  same  graph. 

Table  I  summarizes  the  longitudinal  velocities  at  room  temperature  for  the 
various  plastics  discussed  in  this  section. 


Table  1 

LONGITUDINAL  VELOCITY  IN  COMMON  PLASTICS  AT  25 °C 


Material 

Density 

g/cc 

Frequency  range 

Me 

Velocity 
mm/ usee 

Poly-Penco  nylon 

1.15 

1.5-10.0 

2.72 

Delrin  Acetal  nylon 

1.43 

1. 5-3.0 

2.48 

Polyethylene 

low  density 

0.916 

1.5-4.33 

1.99 

high  density 

0.964 

1.5-10.0 

2.55 

Plexiglass 

1.19 

1.5-10.0 

2.73 

Teflon 

2.19 

1.5-4.33 

1.40 
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2.  Longitudinal  Velocity  in  Composite  Nose-Cone  Materials 

Another  group  of  materials  studied  consists  of  typical  nose-cone  materials 
which  are  generally  either  composites  or  laminates.  In  some  of  these  reentry 
vehicle  materials,  the  measurement  of  ultrasonic  velocities  is  very  difficult 
due  to  excessive  attenuation  and  distorcion  of  the  signal  and  the  results  ob¬ 
tained  should  be  viewed  with  caution  as  far  as  absolute  velocity  values  are 
concerned.  In  the  graphs  that  follow,  most  of  the  data  are  fairly  complete  and 
self-consistent,  but  the  cases  where  the  results  are  suspect  will  be  pointed 
out.  These  situations  usually  result  from  the  fact  t’  in  many  of  the  com¬ 
posite  materials  sampling  errors  arise  because  of  significant  luaterial  inhomo¬ 
geneities  and  nonuniformities.  In  some  of  the  layered  materials  it  is  often 
difficult  to  assign  a  particular  direction  to  the  lamination  layers,  or  to 
reproduce  a  particular  orientation,  so  that  sampling  errors  again  become  impor¬ 
tant. 

The  material  represented  in  figure  11  is  Castable  124  (C  124) ,  a  light- 
yellow,  clear  material  with  a  density  of  1.23  g/cc.  As  can  be  seen,  the  data 
are  insensitive  to  the  different  frequencies  utilized  and  are  fairly  consistent 
so  that  appreciable  dispersion  is  not  present  in  this  material  for  the  frequency 
range  of  1.5  -  3.0  Me.  The  velocity  curve  for  this  material  was  not  corrected 
for  thermal  expansion  because  of  the  difficult;  in  obtaining  reliable  expansion 
data.  When  exposed  to  temperatures  higher  than  75°C,  the  specimen  showed 
significant  dimensional  instabilities  when  remeasured  at  room  temperature. 

These  dimensional  changes,  which  were  as  high  as  3  percent  for  exposure  to 
225*C,  more  than  offset  any  error  resulting  from  not  applying  a  thermal  expan¬ 
sion  correction. 

The  next  two  graphs  show  the  velocities  measured  in  two  similar  materials, 
Avcost  I  (yellow,  density  1.10  g/cc)  and  Avcoat  19  (green,  density  1.07  g/cc). 
Attenuation  in  these  samples  Increased  greatly  at  elevated  temperatures  so  that 
measurements  could  be  made  only  by  using  thinner  samples  and/or  lower  frequen¬ 
cies.  The  velocity  curve  in  Avcoat  19  (Fig.  12)  was  checked  by  three  different 
techniques.  Measurements  were  first  made  at  room  temperature  by  measuring 
transit  times  through  two  different  thicknesses  of  material  (comparison  method), 
but  the  method  did  not  work  well  at  higher  temperatures  because  of  excessive 
attenuation  through  the  thicker  (19  mi)  sample-  Therefore,  the  echo  technique 
deecrlbed  in  Section  II  waa  used  by  placing  the  thinner  sample  (3.5  mm)  between 
two  different  aluminum  transmission  blocks.  This  method  gave  good  agreement 
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Figure  11.  Longitudinal  Wave  Velocity  ve  Temperature  In  Caetable  124 

(denelty  1.23  g/cc) 
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Figure  12.  Longitudinal  Wave  Velocity  v*  Teaperature  in  Avcoat  19 

(density  1.07  g/cc) 
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with  the  values  previously  determined  at  room  temperature  for  Avcoat  19  and  was 
therefore  used  for  the  measurements  at  50  and  75 °C.  To  check  the  velocity 
values  obtained  by  the  echo  technique,  and  to  extend  the  temperature  measure¬ 
ments,  the  single  sample  method  was  used  to  obtain  room  and  high  temperature 
values.  At  high  temperatures  the  specimen  became  soft  and  somewhat  distorted 
even  under  the  approximately  1/3  psi  pressure  from  the  bonding  clamp.  However, 
dimension  measurements  on  the  sample  at  room  temperature  after  exposure  to  160°C 
showed  a  thickness  decrease  of  only  0.4%,  so  that  the  high  temperature  values 
should  be  fairly  insensitive  to  dimensional  instabilities.  The  faces  of  the 
sample  which  were  not  protected  by  the  buffer  rods  darkened  significantly  at 
160°C,  but  the  effect  was  only  observed  at  the  edges.  A  signal  could  not  be 
reliably  detected  at  higher  temperatures  in  this  material,  so  the  velocity  data 
are  limited  to  this  tempera  cure. 

Avcoat  I  (Fig.  13)  exhibited  behavior  similar  to  that  of  Avcoat  19  at  higher 
temperature  in  that  the  material  became  somewhat  soft,  and  the  signal  was  great¬ 
ly  attenuated.  However,  measurements  were  possible  to  250°C  using  the  single 
sample  method,  which  checked  reasonably  well  with  the  extrapolation  of  the 
values  obtained  by  the  comparison  method  at  lower  temperatures. 

After  exposure  to  250°C,  the  sides  of  the  Avcoat  I  specimen  which  were  not 
protected  by  the  oil  film  and  buffer  became  excessively  charred.  However  this 
effect  was  observed  only  for  the  exposed  sides;  the  center  of  the  specimen 
appeared  to  be  unaffected  except  for  some  darkening.  This  material  also  showed 
greater  dimensional  change  than  did  Avcoat  19-  After  exposure  to  250“C,  the 
thickness  of  the  specimen  was  2  percent  lower  than  the  original  value. 

Figure  14  shows  the  temperature  dependence  of  velocity  in  chopped  nylon 
phenolic  (CMP),  density  1.68  g/cc.  The  axis  of  propagation  was  approximately 
perpendicular  to  the  reinforcing  material  (c-direction) .  As  indicated  in  the 
graph,  dispersion  is  not  apparent  in  this  frequency  range,  although  the  measure¬ 
ments  are  somewhat  limited  because  the  attenurtion  increased  rapidly  with  fre¬ 
quency. 

Thr*  velocity  through  pyrolytic  graphite  (density  2.19  g/cc)  was  measured 
in  the  a  and  c  directions  (parallel  and  perpendicular  to  the  reinforcing  layers , 
respectively),  using  the  single  sample  method.  As  shown  in  figure  15,  there 
appears  to  be  no  velocity  dependence  upon  frequency;  however,  it  is  noted  that 
there  is  about  a  2  percent  variation  among  the  data  points  due  to  the  small 
thickness  of  the  sample  available.  The  material  contained  alternate  thin 
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Figure  13.  Longitudinal  Wave  Velocity  vs  Temperature  in  Avcoat 

(density  1.10  g/cc) 
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Figure  15.  Longitudinal  Wave  Velocity  vs  Teaperature  in 
pyrolytic  Graphite  (density  2.19  f/cc) 
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(0.05  -  0.10  mm)  and  thick  (0.2-0. 6  mn)  layers  so  that  the  highest  frequency 
used  here  corresponded  to  a  minimum  wavelength  slightly  greater  than  the  maxi¬ 
mum  layer  thickness.  For  values  of  wavelength  much  less  than  the  layer  spacing 
(corresponding  tc  frequencies  of  5.0  and  10.0  Me)  the  output  waveform  was  dis¬ 
torted  to  the  point  that  reliable  velocity  information  could  not  be  obtained. 

The  curves  shown  in  figure  15  are  corrected  for  thermal  expansion  of  the  sample 
using  average  values  of  the  linear  coefficients  of  thermal  expansion  for  the 
two  orientations. 

Preliminary  measurements  on  several  laminated  materials  definitely  estab¬ 
lished  the  conclusion  that  the  velocity  of  propagation  is  a  sensitive  function 
of  the  angle  of  propagation  through  the  layers.  Three  materials  were  conse¬ 
quently  studied  in  detail  with  the  specific  aim  of  determining  the  angular  de¬ 
pendence  of  velocity.  For  these  measurements  the  single  sample  technique  was 
used  rather  than  the  comparison  method  because  of  the  difficulty  in  obtaining 
two  samples  with  the  same  lamination  angle.  Wave  velocities  in  these  materials 
were  studied  as  a  function  of  lamination  angle  at  several  temperatures  between 
25°  and  125°C. 

The  range  of  frequencies  previously  considered  was  chosen  on  the  basis  that 
the  period  of  the  waves  corresponded  approximately  to  pulse  lengths  of  interest. 
The  variation  of  frequency  with  layer  thickness  discussed  above  indicated 
negligible  frequency  dependence  of  velocity,  within  experimental  variation,  for 
the  materials  .studied  and  the  frequency  range  covered.  For  practically  all  the 
above  materials  the  wavelength  was  much  larger  than  the  layer  spacing  at  3.0  Me. 
This  situation  is  trua  for  the  following  three  materials  also,  and  all  measure¬ 
ments  are  reported  for  a  frequency  of  3.0  Me. 

Figures  16  and  17  show  the  angular  dependence  of  the  longitudinal  velocity 
in  General  Electric  Phenolic  Fiber  Glass  (PFG,  density  of  1.91  g/cc).  The 
majority  of  PFG  angular  measurements  were  made  at  30*C,  although  measurements 
at  angles  of  1#,  23*,  44*,  and  89*  were  also  made  at  temperatures  of  50*C,  75*C, 
and  75*C.  These  measurements  show  the  consistent  shape  of  the  velocity  vs  lami¬ 
nation  angle  curve  as  the  temperature  is  varied. 

The  change  in  the  velocity  as  a  function  of  lamination  angle  in  Phenol 
Carbon  (PC,  density  1.48  g/cc)  shown  in  figures  18  and  19  is  less  than  the 
change  observed  in  PFG.  For  this  reason  there  is  more  uncertainty  in  drawing  a 
unique  curve  to  fit  the  data  points.  The  chosen  curve,  howeve-  appears  to  fit 
all  data  points  better  than  any  ocher  simple  curve. 
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Figur*  17.  Loo|ltudlaal  Wava  Velocity  vs  Tsapsratura  with  Lamination 
Angls  as  Paramatar  la  G8  Phaaolic  Plbsr  Glass 
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Figure  18.  Longitudinal  Wave  Velocity  vs  Lamination  Angles  with  Temperature  as  Parameter 

in  Phenolic  Carbon  (density  1.48  g/cc) 
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In  Tape  Wound  Nylon  Phenolic  (TWNP,  density  1.21  g/cc)  there  is  almost  no 
dependence  of  wave  velocity  on  the  lamination  angle  (figures  20  and  21).  In 
this  case  a  straight  line  is  used  to  fit  the  data  points  in  figure  20.  However, 
a  curve  with  the  same  shape  as  the  PFG  and  the  PC  curves,  only  with  less  verti¬ 
cal  variation,  fits  the  TWNP  data  equally  well. 

The  similar  shapes  of  the  PFG,  the  PC,  and  the  TWNP  velocity  vs  lamination 
angle  curves  led  to  an  attempt  to  predict  these  curves  theroretically.  Horio 
(Ref.  7)  discusses  the  elastic  moduli  of  layered  paper  from  general  considera¬ 
tions  of  elasticity.  With  hi^  approach  he  derived  an  expression  giving  the 
stress  to  st~ain  ratio  (Young's  modulus)  for  an  arbitrary  direction  with  respect 
to  the  plane  layers.  He  found  that  the  expression  thus  obtained  fitted  his 
experimental  data  quite  well. 

Since  the  wavelength  of  the  ultrasonic  waves  in  the  present  case  is  gener¬ 
ally  much  larger  than  the  layer  spacing,  his  approach  should  be  applicable  to  the 
present  experimental  data.  At  3.0  Me  the  number  of  complete  layers/mm  is  about 
4.5  for  PFG,  3.0  for  PC  and  2.3  for  TWPN.  Here,  a  complete  layer  is  considered 
to  be  a  layer  of  the  fiber  glass  cloth  and  a  layer  of  phenolic.  The  wavelength 
(defined  as  1  «  where  c  is  the  velocity  and  f  is  the  frequency)  at  normal 
incidence  and  at  3.0  Me  for  each  case  is  1.1,  1.25,  and  0.9  mm,  respectively. 
Following  Horio' s  approach,  the  longitudinal  wave  velocity  vs  lamination  angle 
curves  are  given  by 
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where  0  is  the  angle  between  the  normal  and  the  plane  of  the  layers,  and  cq  and 
ego  are  the  velocities  for  parallel  and  normal  angles  of  incidence,  respectively. 
With  equation  (2)  and  a  knowledge  of  the  wave  velocity  at  0  -  0*  and  90°,  the 
velocity  for  any  angle  of  incidence  can  be  estimated  to  within  2  percent.  The 
curve  thus  obtained  is  shown  for  comparison  with  the  actual  curve  of  PFG  in 
figure  22.  The  results  of  this  work  indicate  that  equation  (2)  can  be  applied 
to  any  material,  such  as  pyrolytic  graphite,  if  the  wavelength  is  very  much 
larger  than  the  layer  spacing,  and  if  the  velocities  at  angles  of  0*  and  90* 
with  respect  to  the  lamination  layers  are  known. 
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Figure  20.  LoagiCudlaal  Wave  Velocity  va  Lamination  Angle  with  Temperature 
aa  Parameter  in  Tape  Wound  Phenolic  Nylon  (density  1.21  g/cc) 
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Figure  21 ,  Longitudinal  Wave  Velocity  va  Temperature  with  Lamination  Angle 
ae  Parana ter  in  Tape  Wound  Phenolic  Nylon 
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DIFFERENCE 


LAMINATION  ANGLE,  dtgrtu 

Figure  22.  Looeituuinel  Wave  Velocity  in  Phenolic  Fiber  Glee* 
and  the  Difference  Between  the  Observed  end  Predicted  Values 
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Less  complete  velocity  data  were  also  obtained  for  a  few  other  layered  nose- 
cose  materials.  Rather  limited  information  was  obtained  for  two  different  con¬ 
figurations  of  AVCO  Phenolic  Fiber  Glass  (AVCO  PFG).*  When  the  angle  of  incidence 
of  the  wave-5  is  perpendicular  to  the  reinforcing  layers  of  glass  (Fig.  23,  lc  r 
curve)  it  'as  found  that,  for  frequencies  higher  than  about  1  Me,  the  material 
apparently  acts  as  a  low-pass  filter.  When  freq  encies  higher  than  this  are 
transmitted  into  the  sample,  the  received  signal  is  observed  to  have  a  lower 
frequency  than  that  of  the  input.  At  a  frequency  of  1  Me  or  below,  the  trans¬ 
mitted  and  received  signals  are  identical  in  frequency.  However,  the  cutoff 
frequency  of  1.0  Me  is  not  sharp,  since  the  output  frequency  for  3.0  Me  input 
varies  from  1.2  to  2.2  Me,  depending  on  the  particular  sample.  The  room  tempera¬ 
ture  points  shown  for  the  lower  curve  (c-direction)  in  figure  23  were  obtained 
for  several  different  samples.  It  was  found  that  even  though  the  output  fre¬ 
quencies  at  3.0  and  4.0  Me  are  significantly  lower  than  the  input,  the  velocity 
computed  by  the  single  sataple  technique  agreed  to  within  1  percent  of  that  com¬ 
puted  by  the  comparison  method  when  the  first  peaks  were  used  as  references.  For 
this  reason  the  temperature  measurements  were  extended  by  the  single  sample 
method,  because  of  excessive  attenuation  of  the  acoustic  signals  at  high  tempera¬ 
ture. 

The  top  curve  shown  in  figure  23  represents  the  velocity  in  the  a-direction 
through  AVCO  PFG.  As  for  the  lower  curve,  room  temperature  points  were  obtained 
by  using  several  different  samples,  and  the  temperature  measurements  were 
e:  .ended  by  using  one  thin  sample  For  this  orientation,  the  output  frequency 
was  slightly  lower  than  the  input,  ranging  from  2.0  to  2.5  Me  for  3.0  Me  input 
frequencies.  For  both  orientations  the  reduced  frequency  efect  appeared  to  be 
geometric,  since  the  thinner  samples  always  gave  higher  output  frequencies  for 
a  given  input  than  did  the  thicker  samples.  For  this  material  the  number  ot 
complete  layera/mm  is  ''-A,  so  that  a  frequency  of  4.0  Me  corresponds  to  a  wave¬ 
length  of  approximately  t  ree  vimes  a  complete  layer  spacing. 

Numerous  efforts  to  take  meaningful  velocity  measurements  in  Oblique  Tape 
Wound  Refrasil  (0TWR)**  were  somewhat  unsuccessful  in  that  very  erratic  results 
were  obtained.  As  in  the  case  of  the  othar  onisotropic  and  inhomogeneous 


•Density  1.71  g/cc 
•♦Density  1.54  g/cc 

3b 


AFWL-TS-65-188 


Figure  2i.  Loogitud.aal  Wave  Velocity  vs  Temperat<  »  in  AVCO  Fiber  Glass 

(density  1.71  g/cc) 


AFWL-TR-6 5-188 


materials,  different  samples  were  cut  so  as  to  study  the  effect  of  the  angle  of 
Incidence  on  the  velocity  of  propagation.  Since  the  layers  in  OTWR  are  very 
irregular  (not  planar),  no  single  direction  could  properly  be  described  as  either 
perpendicular  to,  or  parallel  with,  the  layers.  Figure  24  shows  the  curves 
representing  the  velocities  in  the  perpendicular  and  parallel  directions.  It  is 
interesting  to  note  the  fairly  large  difference  in  the  velocity  between  0.b67 
and  4.0  Me  for  the  upper  curve.  The  situation  is  similar  for  the  bottom  curve, 
with  velocities  for  0.667  and  1.5  Me  waves  being  2.003  and  2.036  mm/ysec, 
respectively.  To  check  the  consistency  of  this  effect,  the  velocity  was  measured 
at  another  orientation  (^45°)  of  OTWR,  and  the  same  effect  was  observed  (the 
velocity  at  1.5  Me  is  2.732  mm/ysec  and  at  3.0  Me  it  is  2.890  mm/usec) .  These 
effects  suggest  that  this  range  of  wavelengths  must  be  comparable  to  the  layer 
spacings.  It  was  found  that  the  spacing  is  1  mm  in  OTWR,  and  in  the  frequency 
range  of  0.667  to  4.0  Me  the  wavelength  ranges  from  5  to  0.9  mm  for  the  a- 
directlon,  and  from  3  to  0.7  mm  for  the  c-direction.  As  the  wavelength  is  varied 
through  this  region,  it  would  be  expected  that,  because  of  geometric  effects, 
the  velocity  should  change  significantly. 

OTWR  also  behaved  as  a  low  pass  filter  with  a  cutoff  frequency  of  approxi¬ 
mately  1  Me  for  the  a-direction  (all  the  reported  frequencies  were  passed  in  the 
c-direction).  The  geometric  effect  of  the  cutoff  frequency  is  predominant  for 
this  material.  It  was  found  that  a  3-mm  piece  of  OTWR  (propagation  in  the  a- 
dlrection)  passe  '  all  frequencies  to  4,0  Me.  For  a  thicker  piece  (6  mm)  the 
output  frequency  consisted  of  a  small  amplitude  3.0-Mc  wave  superimposed  on  a 
higher  amplitude  wave  whose  frequency  was  approximately  0.67  Me.  For  a  still 
thicker  piece  (12mm)  the  3.0-Mc  carrier  was  completely  absent,  and  the  pre¬ 
dominant  frequency  was  0.67  Me.  The  same  effect  occurred  for  6-Mc  and  9-Mc 
»  Inputs  and  !s  thought  to  be  due  to  geometric  effects.  This  effect  has  not  beer. 

observed  in  moat  materials,  except  to  a  lesser  extent  in  AVCO  Phenolic  Fiber 

Glass. 

One  thin  piece  of  RAD-60  (density,  1.41  g/cc)  w*s  available  for  use  as  a 
sample.  A  measurement  of  velocity  through  this  material  was  made  as  a  function 
of  temperature  using  the  single  sample  technique.  Figure  25  shows  the  relation¬ 
ship  between  velocity  and  teeqperature  to  ^200*C.  It  is  interesting  to  note  the 
curvature  of  the  velocity  curve  in  this  region.  For  the  materials  previously 
studied  it  was  found  that  the  velocity  curve  normally  showed  either  xero  or 
positive  curvature  near  200*C  whereas  the  curvature  for  this  specimen  was 
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Figure  24.  Longitudinal  W«ve  Velocity  vs  Tenpers '”-e  in  Oblique 
Tape  Wound  tofresil  (density  1.54  g 
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25.  Longitudinal  Wave  Velocity  vs  Temperature  in  Rad  60 
(density  1.41  g/cc) 
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negative  in  the  same  region.  This  material  also  exhibited  some  thermal  degrada¬ 
tion  or  physical  change  when  exposed  to  200°C,  since  a  redetermination  of  the 
velocity  at  room  temperature  showed  a  decrease  of  ^2  percent. 

Figure  26  shows  the  velocity  as  a  function  of  temperature  in  Rad  58E  (densi¬ 
ty  1.26  g/cc),  which  is  a  composite  of  85  percent  silica  fibers  and  15  percent 
phenolic  resin.  The  attenuation  in  this  material  is  fairly  high,  and  for  that 
reason  the  curve  is  limited  to  approximately  135°l. 

The  ultrasonic  velocities  in  the  nose-cone  materials  discussed  in  this 
section  are  given  in  Table  II.  The  values  all  correspond  to  25  °C,  and  the 
frequency  range  and  direction  of  propagation  with  respect  to  the  layers  are 
indicated  where  applicable. 

Fused  Silica  (density  1.93  g/cc)  is  also  listed  under  this  group  of  materi¬ 
als  although  it  rightfully  belongs  in  a  class  of  Its  own.  As  in  most  of  the 
materials  previously  discussed,  the  velocity- temperature  variation  in  this 
materiel  was  obtained  by  using  an  oil  bath  for  temperature  regulation.  It 
was  observed,  however,  that  the  fused  silica  used  in  these  measurements  readily 
absorbed  the  oil  from  the  temperature  bath.  Since  the  effect  of  this  would  be 
to  increase  the  overall  density  of  the  specimen,  one  would  expect  a  si  ght  shift 
in  the  velocity.  This  shift  would  occur  because  the  effect  of  tue  oil  would  be 
to  slightly  perturbate  the  stress-strain  matrix  and,  consequently,  to  modify  the 
equations  of  motion.  However,  measurements  on  the  longitudinal  velocity  both 
before  and  after  the  •sample  had  been  immersed  in,oil  showed  no  observed  devia- 
t '.on  due  to  oil  absorption,  and  the  temperature  variation  of  the  longitudinal 
velocity  thus  obtained  is  shown  in  figure  27.  The  data  shown  in  figure  27  are 
corrected  for  thermal  expansion  and  indicate  that  the  velocity  Increases  slight¬ 
ly  with  temperature  to  at  least  125*C,  contrary  to  the  behavior  in  most  materials. 

before  concluding  this  section,  some  discussion  concerning  the  accuracy  of 
the  measurements  is  in  order.  The  sources  of  error  inherent  in  this  method  of 
measuring  ultrasonic  velocities  fall  into  four  major  categories.  These  are: 

(1)  errors  in  sample  thickness  measurements,  (2)  transit  time  measurement  errors, 
(3)  errors  in  temperature  measurements,  and  (4)  systematic  errors. 

One  problem  of  sample  preparation  was  in  machine  shop  technique.  The  early 
samples  used  in  this  study  were  cut  on  a  lathe.  This  procedure  was  found  satis¬ 
factory  for  the  thicker  samples  and  typically  gave  samples  varying  in  thickness 
by  no  more  than  0.05  percent.  With  the  thinner  and  mors  flexible  samples,  how¬ 
ever,  the  thickness  variation  was  observed  to  be  as  high  ar  +  1.7  percent  in  some 
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Table  II 

LOGITUDINAL  VELOCITY  IN  NOSE-CONE  MATERIALS  AT  25°C 


Density 

Frequency  Range 

Velocity 

Material 

r/cc 

Me 

mm/usec 

Castable  124 

1.23 

1. 5-4.0 

2.45 

Avcoat  I 

1  10 

1. 5-3.0 

2.46 

Avcoat  19 

1.07 

1. 5-4.0 

2.14 

Chopped  Nylon  Phenolic 
c-direction 

1.20 

1. 5-3.0 

2.55 

Pyrolytic  Graphite 
a-direction 

2.19 

1. 5-4.0 

4.81 

c-direction 

1. 5-4.0 

3.42 

G.  E.  Phenolic  Fiber  Glass 
a-direction 

1.91 

3.0 

4.34 

c-direction 

3.0 

3.33 

Phenolic  Carbon 
a-direction 

1.48 

3.0 

4.10 

c-direction 

3.0 

3.82 

Tape  Wound  Nylon  Phenolic 
a-direction 

1.21 

3.0 

2.83 

c-direction 

3.0 

2.77 

AVCO  Phenolic  Fiber  Glass 
a-direction 

1.71 

3. 0-4.0 

4.35 

c-direction 

3. 0-4.0 

2.72 

Oblique  Tape  Wound  Refrasil 
a-direction 

1.54 

3. 0-4.0 

3.36 

c-direction 

3. 0-4.0 

2.93 

Rad  60 

1.41 

3. 0-4.0 

2.97 

Rad  58B 

1.26 

3.0 

2.70 

Fused  Silica 

1.93 

1. 5-3.0 

4.42 
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Figure  27.  Longitudinal  Wave  Velocity  vs  Temperature  in  Fused  Silica 

(density  1.93  g/cc) 
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Another  problem  associated  with  sample  thickness  was  that  so  ^e  of  the 
plastics  studied  in  this  work  showed  dimensional  instability.  Teflon,  in  par¬ 
ticular,  exhibited  this  characteristic.  Teflon  samples  initially  received  from 
the  machine  shop  were  very  uniform  in  thickness,  but  after  repeated  c,  '.linj.  to 
high  temperatures  it  was  noted  that  the  velocity  values  were  becoming  increas¬ 
ingly  erratic.  When  the  samples  were  remeasured,  it  was  discovered  that  they 
had  become  noticeably  thicker  and  that' opposite  faces  were  no  longer  parallel. 
This  phenomenum  was  at  least  partially  due  to  the  relaxation  of  strains  ->osed 
during  the  machining  process.  After  this  discovery  uniform  samplts  ..ere  tain- 
ed  by  using  a  milling  process.  Additionally,  frequent  checks  we: °  mad^  of  t  e 
sample  dimensions  between  temperature  measurements  to  reduce  the  error  re?ul 
ing  from  this  cause. 

The  computed  velocities  for  t  2  common  plastics  discussed  in  the  . irst  s 
tion  were  corrected  for  thermal  expa  ;sion  using  handbook  values  for  the  co¬ 
efficients  of  linear  expansion.  In  chose  cases  in  which  the  handb<  jk  ga\e  a 
range  of  values  the  median  value  was  used  as  indicated  in  Table  III.  Assi  ning 
that  the  correct  values  for  the  samples  used  were  within  the  published  rat  ge, 
the  maximum  possible  velocity  error  due  to  the  use  of  median  values  would  be 
+  0.5  percent  for  polyethylene  at  100®C.  This  possible  source  of  error  would 
of  course  become  less  significant  as  the  temperature  was  lowere',  becoming  non¬ 
existent  at  room  temperature. 

When  thermal  expansion  data  were  available  for  the  nose-cone  materials,  they 
were  used  to  correct  the  velocity  temperature  curves.  In  those  case?  where  this 
information  waa  not  available,  the  thermal  axpanaion  was  obtained  by  measuring 
the  length  change  as  a  function  of  temperature.  However,  in  comparing  rhe 
values  computed  in  this  way  with  several  available  handbook  values,  it  was  found 
that  tha  agreement  in  the  velocity  curves  was  within  0.5  parcant  at  1C0*C, 
conai  ering  tha  difference  between  the  reported  and  the  experimental  valuer  of 
thermal  axpanaion.  fable  III  lists  the  valuea  of  thermal  expansion  coefficients 
ue ad  for  the  previously  mentioned  materials.  The  letters  a  and  c  after  the 
number  indicate  the  a  and  c  directions  respectively. 

Errors  in  time  measurements  might  srlst  from  three  sources.  Two  of  thest, 
srrors  from  tht  tims  mark  gensrstor  and  tha  oecilloscop may  well  be  considered 
negligible  in  coaq>*rlson  with  other  sources  of  error.  Tht  most  important  source 
of  tims  error  is  due  to  the  interpolation  method  of  measuring  time  intervals  and 
arises  because  of  operator  technique.  The  oscilloscope  dslay-potenliomscer  was 
reed  to  three  decima'  places  by  visual  interpolation  batveen  the  finest 
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Table  III 

COEFFICIENTS  OF  LINEAR  EXPANSION  FOR  SEVERAL 
COMMON  PLASTICS  AND  NOSE-CONE  MATERIALS 


Material 

Nylon  (Poly  Penco  and  Delrin  Acetal) 

Polyethi^ne  (high  and  low  density) 

Plexiglas 

Teflon 

Avcoat  I 

Avcoat  19 

Chopped  Nylon  Phenolic 
F, roly tic  Graphite 

Phenolic  Fiber  Glaaa 

Tape  Wound  Nylon  Phenolic 

AVCO  riber  Glass 

Oblique  Tape  Wound  Ee-fraail 

Rad  60 
Rad  58B 
Fused  Silica 


Coefficients  of  Expansion 
(10  6mo/nnn  *C) 


90 

55 

90 

150 

63  (a) 

18  (a) 
26  (c) 

13  (a) 
45  (c) 

90  (a) 
130  (c) 

20  (a) 

13  (c) 

10  (a) 
20  (c) 

14 
40 


1.1 


* 
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gradations  on  the  scale.  From  reproducibility  measurements  between  operators, 
it  is  doubted  that  the  error  in  making  any  individual  time  measurements  ex¬ 
ceeded  0.01  psec.  Since  the  difference  between  elapsed  times  through  two 
samples  was  typically  in  the  range  from  4  to  6  usee,  the  relative  uncertainty 
in  time  measurements  was  then  not  more  than  0.5  percent. 

The  third  type  of  error  enumerated  above  concerns  temperature  measurements. 
Temperatures  were  measured  with  either  a  thermometer  or  a  thermocouple  and  were 
recorded  to  the  nearest  degree.  The  instrument  used  was  calibrated  to  the  near¬ 
est  degree  over  the  range  of  interest.  The  typical  effect  of  a  temperature 
error  can  easily  be  determined  from  the  linear  graph  of  figure  8,  where  the  func¬ 
tional  relationship  between  velocity  and  temperature  is  seen  to  be  approximately, 

c  -  2.58  -  0.0078(T-20)  (3) 

where  c  is  the  longitudinal  velocity  in  mm/usec,  and  T  is  in  °C.  An  error  of 
one  half  degree  in  temperature  would  cause  an  apparent  error  of  0.004  mm/usec, 
which  is  a  relative  error  of  approximately  0.2  percent  at  the  high  temperature 
end  of  the  curve. 

The  remaining  sources  of  error  other  than  those  discussed  above  constitute 
the  fourth  category.  The  errors  in  this  group  are  systematic,  and,  consequent¬ 
ly,  a  numerical  value  cannot  readily  be  assigned.  The  major  premise  of  this 
experimental  technique  is  thu.  the  acoustic  units  can  be  assembled  accurately 
enough  so  that  units  assembled  using  different  samples  of  the  same  material  will 
differ  ultrasonically  only  because  of  the  differences  in  sample  thickness. 

Because  of  the  time  required  to  statistically  determine  this  factor,  and  because 
of  the  accuracy  it  would  require,  it  was  not  feasible  to  determine  this  factor 
exactly.  However,  comparison  with  standards  and  the  results  obtained  on  the 
materials  under  study  Indicate  that  this  error  is  insignificant  in  view  of  the 
limited  accuracy  required. 

Another  possible  source  of  error  was  that  involved  in  choosing  the  refer¬ 
ence  points  on  the  transmitted  snd  received  pulses  between  which  to  measure  the 
elapsed  transit  tims.  The  received  signals  did  not  coincide  exactly  in  shape 
with  the  rransal-ted  signals,  and  some  small  amount  of  judgment  waa  necessary 
to  pick  the  proper  reference  points.  Usual  however,  the  oscillations  were 
so  clear  and  well  defined  that  there  was  no  serious  doubt  about  the  proper 
choice. 
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When  the  velocities  were  computed  by  using  the  transit  times  through  three 
different  samples  of  the  same  composition,  three  interrelated  velocity  values 
were  obtained  for  each  determination.  This  technique  contributed  greatly  to  the 
elimination  of  error  and  gave  assurance  that  reasonably  accurate  values  were 
obtained.  However  when  reasonable  care  was  exercised  in  assembling  the  acoustic 
unit,  the  accuracy  of  the  single  sample  method  was  found  to  agree  quite  well 
with  that  obtained  by  the  comparison  method. 

Because  not  enough  measurements  were  made  of  the  same  quantity  to  permit  a 
statistical  error  analysis,  a  series  of  maximum  error  calculations  was  made. 
These  showed  that  the  maximum  relative  uncertainty  in  velocity  could  be  as  high 
as  4  percent.  Examination  of  the  data  points  for  which  multiple  velocity  deter¬ 
minations  were  made  on  homogeneous  materials  reveals  a  maximum  spread  of  less 
than  3  percent;  it  should  be  recalled,  however,  that  these  points  include  deter¬ 
minations  made  by  different  timing  methods  and  at  different  frequencies. 

Another  factor  to  be  considered  is  that  many  of  the  early  determinations  were 
made  with  samples  which  were  later  found  to  be  irregular  in  thicknees.  When 
homogeneous,  isotropic  samples  of  sufficient  thicknesses  were  measured  it  was 
observed  that  the  spread  in  Velocity  values  rarely  exceeded  1  percent,  for  a 
material  *  assured  at  a  particular  temperature  and  frequency. 

3.  Longitudinal  Velocity  in  Solid  Epoxv  Foams 

This  section  presents  data  concerning  the  longitudinal  bulk  velocity  in 
solid  epoxy  foams  as  functions  of  temperature  and  density.  The  different  densi¬ 
ties  in  this  rigid  material  were  obtained  by  varying  the  amount  of  blowing  agent 
and  epoxy  resin,  thereby  obtaining  the  appropriate  ratio  of  void  to  epoxy  in  :he 
matrix  of  the  resulting  structure  (Ref.  8).  The  resulting  material  exhibited  a 
fine-cell,  uniform,  nondirectional  structure  completely  free  from  grain  and 
characterised  by  excellent  machinability  and  dimensional  stability. 

Since  the  attenuation  of  longitudinal  bulk  waves  in  this  material  greatly 
exceeded  that  of  the  other  materials  studied,  the  acoustic  experiment  was 
slightly  modified.  The  electronic  apparatus  necessary  for  the  generation  and 
detection  of  the  ultraaonically  pulsed  signals  was  similar  to  that  described 
before  with  the  exception  that  an  astplifier  of  about  ?0-db  gain  was  inserted 
between  the  receiving  attenuator  and  the  oscilloscope.  The  tra.  sit  time  through 
the  malarial  was  then  obtained  by  the  single  sample  technique  using  samples 
approximately  5mm  thick.  In  piece  of  the  transformer  oil,  silicone  grease  was 
used  for  the  coupling  agent  because  it  is  more  viscous  and  does  not  diffuse 
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away  from  the  surface  into  the  foam  when  the  sample  is  exposed  to  high  tempera¬ 
tures.  The  entire  assembly  as  used  for  longitudinal  measurements  was  held  to¬ 
gether  by  light  spring  tension  (about  1/3  psi) . 

I 

The  longitudinal  bulk  velocity  is  plotted  as  a  function  of  temperature  with 
density  as  parameter  in  figures  28  and  29.  Thermal  expansion  data  was  not  avall- 
ifble  for  this  particular  foam,  so  the  curves  are  not  corrected  for  thermal  ex¬ 
pansion.  However,  comparison  with  solid  phenolic  (Ref.  9)  and  silicone  (Ref.  10) 
resins  in  the  same  density  range  indicates  that  the  resulting  error  due  to  ex¬ 
pansion  is  less  than  1  percent  at  temperatures  to  150°C#  For  the  densities 
shown  in  figures  28  and  29,  the  samples  had  not  been  previously  heat-treated, 
and  the  same  general  shape  of  the  velocity  curve  is  apparent  for  all  densities. 

Of  particular  interest  is  the  existence  of  an  inflection  point  in  the 
velocity  curve  at  a  temperature  of  about  110°C,  where  there  is  a  definite  change 
in  slope.  Velocity  inflection  points  are  generally  associated  with  some  sort 
of  phase  transition  relating  to  a  change  in  one  of  the  moduli.  Similar  behavior 
was  previously  noted  for  Teflon  and  has  been  observed  in  glass  transition  studies 
of  polyethylene  (Ref.  11).  In  this  case,  a  permanent  change  in  the  color,  den¬ 
sity,  and  velocity  was  observed  when  the  samples  were  exposed  to  temperatures 
higher  than  about  110°C.  To  this  temperature,  the  velocity  and  density  values 
were  reproducible  within  experimental  error.  However,  when  the  temperature 
increased  beyond  110°C,  a  permanent  dimensional  change  occured  which  caused  a 
velocity  decrease  of  5  to  10  percent  and  a  density  increase  of  5  to  15  percent 
at  room  temperature.  The  actual  values  were  a  function  of  the  density  and  v 
the  temperature  to  which  the  samples  had  been  exposed,  the  lower-density  samples 
showing  greater  dimensional  changes.  The  color  also  changed  from  the  original 
neutral  to  various  shades  of  brown,  depending  on  the  temperature  to  which  the 
sample  had  been  exposed. 

Figure  30  shows  the  bulk  velocity  variation  with  temperature  of  two  high 
density  foams  that  had  been  previously  heat  treated  and  loaded  with  different 
amounts  of  aluminum.  The  aluminum  filler  used  for  loading  was  finely  ground  to 
an  average  particle  size  of  less  than  25  microns,  thus  insuring  the  nondirection- 
al  quality  of  the  foam.  The  bottom  curve  corresponds  to  a  light  gray  epoxy 
foam  (density  0.371  g/cc)  prepared  in  the  same  fashion  as  the  lower  density 
unloaded  foams  but  filled  with  8  percent  aluminum  by  weight.  The  top  curve 
represents  a  darker  gray  material  (density  1.068  g/cc)  prepared  by  mixing  equal 
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P^rts  by  weight  of  the  finely  powdered  8  percent  leaded  foam  *"d  solid  epoxy 
resin  (Ref.  8).  Both  materials  had  been  previously  heat  treated  at  250*C  for 
several  hours,  and  it  is  noted  that  an  inflection  point  in  the  velocity  does 
not  exist  for  either  material.  The  velocity  and  density  for  both  materials 
were  reproducible  within  experimental  error  at  room  temperature  after  the 
samples  had  been  repeatedly  cycled  to  high  temperatures.  High  temperature  in¬ 
formation  for  the  lower  density  material  was  limited  to  15Q*C  because  of  ex¬ 
cessive  attenuation  of  the  signal  at  higher  temperatures. 

Figure  31  shows  the  attenuafior  in  db/cm  plotted  as  a  function  of  tempera¬ 
ture  for  the  0,325  g/cc  sample.  Th*  "ttenuation  is  calculated  from  the  relation, 


a 


X2  -  Xj  -  c 

— 


(4) 


where  Xj  and  X2  are  the  attenuator  settings  in  db'.  for  a  standard  amplitude 
signal  with  and  without  the  sample  Inserted,  and  d  is  the  thickness  of  the 
sample.  C  is  a  correction  factor  which  accounts  for  the  impedance  mismatch  at 
the  two  buffer-sample  Interfaces  and,  ss  shown  In  Appendix  V,  is  given  by 


C 


20  log 


4Zj 


7. 


(5) 


Zj  end  Z?  being  the  characteristic  acoustic  impedances  of  the  sample  and  alumi¬ 
num  buff  is.*  Equation  (4)  do«»  not  account  for  Inter fecial  effects  such  ss  the 
Impedance  of  the  silicone  grease  or  the  reproducibility  of  the  bonding  surface. 
The  first  effect  was  found  to  be  negligible  by  comparison  with  attenuation 
measurements  obtained  by  using  two  samplss  of  different  thicknesses.  In  the 
latter  case  the  correction  factor  goes  to  zero  sines  the  seme  number  end  types 
of  Impedance  mismatches  are  present  in  both  assemblies.  The  value  of  d  to  be 
used  in  the  resultant  equation  is  the  difference  in  thicknesses  of  the  two 
samplss. 


*The  longitudinal  acoustic  impedance  of  s  medium  is  defined  ss  Zj  -  pjCj,  where 
Pj  end  c j  correspond  to  the  density  and  longitudinal  velocity  In  the  medium.  A 
similar  relation  holds  for  transverse  waves  when  c  represents  the  transverse 
wavs  velocity 
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Figure  31.  Loagltudlnel  Have  Attenuation  ve  Tenperetur*  ir  Zpoxy  Foes 

(density  0.3?-3  g/ee) 
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For  both  of  the  above  methods  the  accuracy  of  the  attenuation  measurements 
was  found  to  rely  on  the  reproducibility  of  component  assembly.  For  samples 
whose  surfaces  are  not  smooth,  the  amount  of  acoustic  energy  transferred  is  a 
sensitive  function  of  the  contact  area  wetted  by  the  coupling  agent.  As  more 
bond  Is  added,  better  transfer  of  erergy  occurs  since  the  area  of  contact,  es¬ 
pecially  for  the  foam  samples,  is  not  limited  to  the  area  of  a  circle  (if  right 
circular  samples  are  used).  This  occurs  because  the  sample  surfaces  are  fairly 
rough  due  to  the  porous  structure  of  the  foam.  Therefore,  the  amount  of  bond 
and  the  applied  pressure  must  be  reproducible  if  accurate  values  of  attenuation 
are  desired.  For  these  particular  measurements,  the  precision  of  the  attenua¬ 
tion  measurements  was  found  experimentally  to  be  less  than  20  percent.  For  the 
reasons  discussed  above  the  attenuation  values  should  be  considered  relative 
rather  than  absolute. 

Of  particular  interest  in  figure  31  is  the  existence  of  an  inversion  in  the 
attenuation  at  about  110’C.  It  will  be  recalled  that  the  velocity  inflection 
occurred  at  this  tertpercture  end  th«t  a  p*manent  change  in  the  slope  of  veloci¬ 
ty  existed  for  higher  temperatures.  The  situation  is  similar  for  the  attenua¬ 
tion  measurements;  however,  the  effect  seems  to  be  much  more  pronounced  and 
consequently  a  better  indication  of  the  presence  of  the  transition. 

Figure  32  shows  the  variation  of  velocity  and  attenuation  as  a  function  of 
frequency  for  the  0.325  g/cc  sample  at  a  constant  temperature  of  22*C.  First, 
it  is  noted  that  extending  the  frequency  to  4.0  Me  increases  the  velocity  about 
20  percent.  This  is  a  common  chain  t of  viscoelastic  materials  in  that 
increasing  the  frequency  is  equivalent  to  lowering  the  temperature  (Refs.  12, 
13).  The  attenuation  is  also  observed  to  increase  with  increasing  frequency, 
and  in  this  case  an  approximately  linear  relationship  exists  with  frequency  for 
both  the  attenuation  and  the  velocity.  Attenuation  and  velocity  measurements 
ware  not  possible  beyond  4.0  Me  with  the  present  technique  because  the  attenua¬ 
tion  increased  very  rapidly  with  higher  frequencies. 

Figure  33  shows  the  variation  of  attenuation  and  velocity  with  density  for 
the  unloaded  samples  at  frequencies  of  C.667  and  1.5  Me.  The  velocity  is  seen 
to  be  essentially  linear  with  density,  monotonically  increasing  in  the  observed 
range.  This  result  is  to  be  expected  since  the  effective  stiffness  should  in¬ 
crease  with  increasing  density,  and,  in  fact,  similar  affects  have  been  observed 
in  density  studies  on  polyethylene  (Ref.  5).  Over  the  same  density  range  the 
attenuation  is  also  found  to  be  linear,  decreasing  with  increasing  density. 
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Figura  32.  Loaf titud Inal  tftv*  Velocity  and  Attenuation  v*  Frequency 
in  Epoxy  tom  of  Daoolty  0.323  g/cc 
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Figure  33.  Loogltudlfial  Ware  Velocity  n  Density  is  Solid  Epoxy  Foes 
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One  further  word  on  the  accuracy  of  the  velocity  measurements  is  required  at 
this  point.  As  noted  in  figures  28  through  30,  the  difference  in  velocity  among 
some  of  the  curves  is  on  the  order  of  1  to  2  percent.  This  magnitude  of  pre¬ 
cision  vouj.6  be  necessary  to  draw  the  individual  curves  of  figures  28  and  29. 

As  previously  mentioned,  the  present  technique  has  been  found  to  give  absolute 
accuracies  on  this  order  under  favorable  conditions,  such  as  sufficient  sample 
thickness,  homogeneous,  linear,  nondispersive  sample®,  and  tha  use  of  thin  coup¬ 
ling  agents  such  as  transformer  oil.  However,  for  these  particular  samples  the 
above  conditions  are  not  entirely  fulfilled  and  the  accuracy  is  consequently 
limited  by  four  major  causes.  First,  as  indicated  previously,  the  thickness  of 
the  interfacial  layer  of  silicone  grease  is  a  critical  factor  in  the  transit 
time  reproducibility  from  one  assembly  to  the  next.  However,  this  error  was 
partially  minimized  by  standardized  techniques  of  assembling.  By  analyzing 
standards  and  using  different  sample  thicknesses,  it  was  found  that  the  pre¬ 
cision  was  well  within  1  percent  when  careful  attention  was  paid  to  the  tech¬ 
niques  of  assembly. 

The  second  error  arises  because  of  the  very  small  thickness  of  the  samples. 
This  can  result  in  a  rather  large  percentage  error  in  the  transit  time  or  can 
magnify  any  error  due  to  sample  nonuiiformity .  Howe’er,  the  former  effect  wos 
small  and  fell  into  the  category  of  the  first  source  of  error  because  of  the 
excellent  machlnablllty  of  these  solid  foams.  It  was  observed  quite  routinely 
thct  the  magnitude  of  the  thickness  variation  was  on  the  order  of  0.25  to  0.50 
percen  for  samples  up  to  5mm  thick.  Thin  tolerance  was  also  found  to  hold 
even  after  the  samples  had  b^en  exposed  to  higher  emperatures  where  dimensional 
changes  might  be  expected  to  occur.  Consequently,  errors  arising  from  the  second 
source  of  error  were  not  deemed  significant. 

The  third  source  of  error  results  from  the  fact  that  relatively  arge  voids 
could  sometimes  be  observed  aloag  the  machined  faces  of  the  prepared  awnple. 

This  is  a  quite  obvious  source  of  error  since  the  velocity  oi  ptcpagatlo  w  uld 
be  a  function  of  the  relative  amount  of  epo*/  encountered  by  the  ultrasonic 
waves.  The  resultant  sampling  error  is  somewhat  indicated  by  the  top  curve  in 
figures  28  and  29.  The  curve,  corresponding  to  the  0.325  g/cc  sample,  consists 
of  points  obtained  from  two  different  samples  cut  from  the  same  block.  The 
measured  denaity  of  each  sample  was  found  to  be  the  same,  and  in  neither  case 
do  the  points  obtained  from  the  two  samples  differ  by  more  then  2  percent. 
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Other  samples  checked  for  sampling  error  yielded  specimens  differing  slightly  in 
density  and  were  consequently  plotted  as  separate  entities  in  figures  28  and  29. 

Hie  last  source  of  error  is  by  far  the  largest  and  depends  upon  the  particu¬ 
lar  definition  used  in  the  calculation  of  velocity.  These  materials  are  obvious¬ 
ly  dispersive,  as  can  be  observed  by  noting  the  velocity  variation  with  frequency 
in  figure  32.  Dispersion  in  elastic  materials  is  generally  due  to  purely  geo¬ 
metrical  consideration  such  as  those  arising  in  thin  rods  or  plates,  whereas  in 
viscoelastic  materials  dispersion  is  produced  by  the  properties  of  the  solid 
(Refs.  14,  1:,,  16).  In  the  latter  case,  high  frequency  waves  travel  faster  than 
low  frequency  ones,  which  is  consistent  with  the  present  experimental  observa¬ 
tions.  But  when  the  elastic  properties  vary  with  frequency,  as  in  a  dispersive 
medium,  the  interpretation  of  the  results  becomes  uncertain,  since  there  is  not 
a  unique  velocity  of  propagation,  and  the  velocity  at  which  energy  is  transferred 
is  the  group  velocity,  c  ,  which  differs  from  the  phase  velocity,  c  ,  by  an 

dcP 

amount  X  -rr—  (X  is  used  here  as  the  wavelength). 

G  A 

As  noted  earlier ,  the  transit  time  required  for  the  pulse  to  traverse  the 
sample  was  computed  by  measuring  from  one  peak  on  the  unrectified  transmitted 
signal  to  the  corresponding  one  on  the  received  signal.  When  the  input  and  out¬ 
put  wave  forms  were  not  distorted  (as  was  the  case  for  most  materials  discussed 
previously) ,  this  method  gave  identical  results  within  the  capabilities  of  the 
technique,  regardless  of  the  peak  used  for  reference.  For  the  foam  materials 
however,  the  calculation  of  the  velocity  was  strongly  dependent  on  which  peak 
was  used  to  measure  the  transit  time.  The  possible  reference  peaks  were  limited 
to  the  first  and  second  positive  (or  negative)  peaks  because  of  subsequent  dis¬ 
tortion  of  the  remaining  peaks  o’  the  output  pulse.  The  velocity  calculated 
usiug  the  second  positive  peak  as  reference  instead  of  the  first  gave  values 
which  were  about  8  percent  lower  at  room  temperature  for  0.667-  and  1.5-Mc  input 
frequencies.  The  reasons  for  this  apparent  inconsistency  is  due  to  the  above 
considerations  and  to  the  fact  that  the  output  frequency  was  in  general  lower 
than  the  input.  For  input  frequencies  of  0.667  to  4.0  Me,  the  output  ranged 
from  0.5  to  2.0  Me.  It  was  also  noted  that  the  output  frequency  slowly  decreased 
with  increasing  temperature  using  these  sane  frequencies,  .'or  an  input  of  0.333 
Me,  the  output  was  comparable,  although  velocity  information  could  not  be  ob¬ 
tained  because  of  extreme  distortion  of  both  input  and  output  waveforms. 

Since  these  foams  are  characteristic  of  visocoelastic  solids,  four  independ¬ 
ent  measurements  would  be  nectssary  to  completely  characterise  the  four  constants 
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of  propagation.  These  constants  consist  of  the  congressional  and  shear  moduli 
•’ssociated  with  elastic  behavior,  plus  the  congressional  and  shear  viscosities. 
To  completely  specify  the  propagation  constants  would  then  require  shear  wave 
velocity  and  attenuation  measurements  as  well  as  the  associated  longitudinal 
measurements.  Shear  wave  measurements  were  attempted  in  the  present  experiments 
by  using  appropriate  coupling  agents  and  extremely  thin  samples.  But  shear  data 
was  not  possible  with  the  present  technique  because  of  the  extreme  attenuation 
of  the  waves.  This  was  to  be  expected  at  these  temperatures  and  frequencies, 
since  the  longitudinal  wave  attenuation  was  high  and  since  shear  waves  are,  in 
general,  attenuated  much  more  than  longitudinal  waves. 

A  Differential  Thermal  Analysis  (DTA)  was  performed  on  the  0.325  g/cc  sample 
(Ref.  17)  and  thermograms  were  obtained  to  approximately  400°C.  It  was  found 
that  a  glass  transition  occurs  at  146°C  with  decomposition  beginning  at  approxi¬ 
mately  200°C  and  complete  decomposition  resulting  at  408°C.  Of  particular 
interest  was  the  existence  of  a  small  curing  exotherm  in  the  region  of  50°  to 
145 °C  which  disappeared  on  repeated  heating.  Hie  material  is  a  thermosetting 
resin  and,  as  such,  crosslinking  of  the  polymer  increases  with  increasing 
temperature.  The  resultant  effect  is  probably  such  that  the  chemical  structure 
changes  with  increasing  cure,  resulting  in  a  permanent  decrease  in  either  the 
shear,  the  compressional  modulus,  or  both,  and  consequently  manifests  itself  in 
a  decrease  of  the  longitudinal  velocity. 

A  word  is  in  order  about  the  behavior  of  the  aluminum  filled  samples. 
Analyses  of  aluminum  have  shown  that  the  longitudinal  velocity,  and  consequent¬ 
ly  the  effective  modulus  for  propagation  of  longitudinal  waves,  s  a  slowly 
decreasing  function  of  temperature  to  at  least  125°C.  One  could  assume  that  the 
properties  of  the  foams  should  become  more  characteristic  or  the  properties  of 
aluminum  as  the  percentage  of  aluminum  increases.  This  assumption  would  in  part 
be  supported  by  the  top  curve  of  figure  30.  However,  as  previously  mentioned, 
the  lower  curve  corresponds  to  a  higher  percentage  of  aluminum  (8  percent  versus 
4  percent  for  the  top  curve)  so  that  the  assumption  is  apparently  invalid  or 
masked  by  more  predominant  characteristics.  If  the  loading  effect  were  compared 
with  the  combined  effects  due  to  density  change,  loading,  and  the  epoxy  itself, 
the  behavior  could  be  definitely  attributed  to  a  particular  cause,  aut  informa¬ 
tion  on  the  variation  of  each  effect  with  the  other  two  held  constant  was  not 
available,  so  that  it  was  not  possible  to  assign  the  differences  in  the  two 
graphs  of  figure  30  to  any  one  effect.  It  seems  reasonable  to  assume  that  the 
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properties  of  epoxy  and  those  inherent  with  the  density  are  responsible  for  the 
difference  in  the  two  graphs,  since  the  lower  curve  (higher  percentage  of  alumi¬ 
num)  corresponds  closely  in  shape  to  those  obtained  from  the  unloaded  samples. 

One  further  word  is  necessary  to  close  the  discussion  on  the  solid  foams. 
Static  tests  on  the  compressive,  flexural,  and  tensile  strengths,  as  well  as 
moduli,  indicate  an  increase  of  approximately  an  order  of  magnitude  in  these 
quantities  when  the  density  of  the  unloaded  samples  was  increased  fourfold 
(Ref.  8).  The  effective  modulus  for  bulk  wave-  propagation  shows  similar  behavi¬ 
or,  increasing  by  almost  8  factor  of  ten  over  the  same  density  range.  As  shown 
in  Appendix  11,  a  relationship  exists  between  the  effective  stiffness  of  a  ma¬ 
terial  and  the  velocity  of  propagation  so  that  comparisons  between  the  static 
and  dynamic  moduli  can  be  made.  The  effective  modulus  is  given  approximately 
by  (Ref.  15,  16), 


A  +  2y  ■  B8  +  4/3y  ■  pc|  (6) 

where  y  and  y  are  the  Lame  constants,  which  relate  stress  to  strain  in  homo¬ 
geneous,  isotropic  materials.  Bs  is  the  adiabatic  bulk  modulus,  p  is  the  densi¬ 
ty,  and  c^  is  the  longitudinal  wave  velocity.  Equation  (6)  is  limited  to  situa¬ 
tions  involving  low  attenuation  and  no  dispersion,  where  the  phase  and  group 
velocities  are  equal.  The  equation  as  presented  here  neglects  the  high  attenua¬ 
tion  in  the  foam,  the  fact  that  the  material  is  dispersive,  and  is  presented 
only  for  comparison  purposes. 

Table  IV  summarizes  the  velocity  through  various  densities  of  the  epoxy 
foams  at  room  temperature. 
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Table  IV 

LONGITUDINAL  VELOCITY  IN  SOLID  EPOXY  FOAM  AT  25°C 


Density 

Frequency 

Logltudlnal  Velocity 

Jt/cc 

Me 

ma/uaec 

0.088 

0.667 

1.08 

1.5 

1.12 

0.154 

0.667 

1.16 

1.5 

?. .  20 

0.166 

0.667 

1.19 

1.5 

1.22 

0.212 

0.667 

1.21 

1.5 

1,25 

0.252 

0,667 

1.28 

1.5 

1.31 

0.261 

0.667 

1.30 

i  5 

1.32 

0.325 

O.t  >7 

1.41 

1.5 

1.45 

3.0* 

1.61 

4.0* 

1.65 

♦Temperature ,  22°C 
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SECTION  IV 
APPLICATIONS 

In  recent  y  ars  many  laboratories  have  concentrated  on  measuring  the  high- 
temperature  properties  of  materials.  This  is  particularly  true  of  space-age 
materials  in  that  they  must  be  capable  of  withstanding  dynamic  vibratory  load¬ 
ings  at  high  temperatures.  Thus,  in  many  practical  applications  it  is  necessary 
to  kn  3  the  values  of  the  elastic  moduli  of  materials  at  high  temperatures. 

These  values  are  necessary  not  only  to  the  design  engineer,  but  also  to  the 
physicist  interested  in'  the  equation  of  state  of  materials  that  make  up  compo¬ 
nents  of  these  systems.  One  method  currently  used  for  making  these  high  tempera¬ 
ture  measurements  is  the  standard  tensile  test  employed  in  the  static  testing  of 
materials  (Ref.  18).  The  specimen  is  enclosed  within  a  furnace  and  connected  to 
the  crossheads  of  a  tensile  machine,  and  stress-strain  plots  are  made  at  elevat¬ 
ed  temperatures.  This  method  of  testing  at  high  temperatures  has  many  short¬ 
comings.  One  is  the  difficulty  in  making  strain  measurements  at  high  tempera¬ 
tures.  In  addition,  the  magnitude  of  the  static  elastic  constant  la  always 
Influenced,  to  some  extent,  by  relaxation  effects  at  elevated  temperatures. 
Superimposed  on  the  purely  elastic  deformation  of  the  specimen  are  deformations 
associated  with  creep  and  elastic  hysteresis.  Thus,  the  values  of  the  elastic 
modulus  obtained  by  static  methods  can  become  very  dependent  upon  the  strain 
rate  used  in  the  test. 

Since  the  propagation  characteristics  of  small-amplitude  elastic  waves  in 
a  solid  are  governed  by  the  elastic  coefficients,  a  measurement  of  the  different 
velcwltles  of  propagation,  such  as  transverse  and  longitudinal  velocities,  allows 
a  convenient  calculation  of  the  elastic  constants.  Ultrasonic  wave  velocity 
measurements  are  easily  made  as  a  function  of  temperature  so  that  the  variation 
of  the  elastic  constants  with  temperature  is  resdily  obtained. 

It.  anisotropic  bodies,  such  as  crystals,  the  number  of  constants  necessary 
to  specify  the  mechanical  properties  of  the  medium  depends  upon  the  particular 
symmetry  system  to  which  the  crystal  belongs.  For  example,  for  cubic  symmetry 
three  coefficients  ere  necessary  to  specify  the  stress  strain  matrix  and  ars 
commonly  designated  as  cu,  cl2.  c44.  In  isotropic  bodies  the  elastic  behavior 
can  be  completely  designated  by  two  elastic  constant*,  X  and  y,  known  as  the 
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Lamd  constants  (Ref.  19).  Xt  Is  shown  In  Appendix  II  that  the  twu  velocities 
of  propagation  In  an  isotropic  can  be  expressed  as 


where  X  and  u  have  been  defined,  c^  is  the "longitudinal  wave  velocity,  cf  is  the 
transverse  velocity,  and  p  is  the  density.  The  Load  constant,  v,  is  identical 
with  the  shear  modulus,  G .  Equation  (7)  applies  to  infinitesimal  displacements 
in  an  infinite  medium,  l.e.,  where  the  wavelength  is  much  smaller  than  the 
longitudinal  or  lateral  dimensions  of  the  specimen.  The  constants  obtained 
from  equation  (7)  are  the  adiabatic  rather  than  the  isothermal  constants,  since 
the  velocity  of  propagation  for  ultrasonic  waves  in  the  low  megacycle  range  is 
under  the  adiabatic  response  of  the  material  (the  thermal  diffusivity  is  Bmall 
at  these  frequencies).  However,  to  compare  with  static  measurements  the  con¬ 
version  to  the  Isothermal  moduli  involves  the  ratio  of  specific  heats  and  is 
usually  a  anal'  correction,  on  the  order  of  a  few  percent  (at  low  pressures). 

Young's  modulus,  E,  the  shear  modulus  G,  the  bulk  modulus  B  and  Poisson's 
ratio  o  can  be  expressed  in  terms  of  the  Land  constants  for  an  isotropic  medium 
as  (Ref.  14) 


E  - 


G  •  U 

B  »  A  +  —y 


a  • 


X 

2(1  +  v) 


(8) 


where  it  ie  understood  thet  ell  of  the  moduli  are  the  adlebetic  moduli.  Subefci- 
tuing  for  the  Lend  constants  in  aquation  (7)  yields  the  following  formulas  in 
Carve  of  the  longitudinal  and  ahesr  wave  velocities 
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in  cgs  or  oka  units. 

Nolle  and  Sleek  (Ref.  4)  discuss  the  mechanical  behavior  of  a  high- loss 
specimen  in  compression  and  shear.  Their  assumption  is  that  for  high-loss  (high 
attenuation  of  acoustic  signals)  the  moduli  must  be  complex,  with  a  real  part 
associated  with  elastic  behavior,  and  an  imaginary  part  which  is  a  function  of 
frequency  and  viscosity.  Following  their  notation,  the  complex  modulus  quanti¬ 
ties  connected  with  the  observed  velocity  and  attenuatlo  are  given  as 


where  p  is  the  density  and  M*  is  the  effective  stiffness  for  a  certain  wave 
motion  whose  velocity  is  given  by  c*.  Direct  solution  of  equatior.3  (10)  gives 


M,  -  p  (c,*  -  C2*) 
-  2pCj  c2 


(ID 


Also,  examination  of  the  usual  complex  exponential  expression  for  propagation 
shows  that  C|  and  c2  are  related  to  the  observed  phase  velocity  c,  to  the  ampli¬ 
tude  constant  a  (db/cm),  end  to  Mj  end  M2  by  the  following  relations 
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1  +  r2 

y 


-2 - 7  c 

i  +  r2  '| 

1  U  +  Y2)2 


M2- 


2Y 


1  - 


(12) 


CIC 

where  y  m  0  1*  the  attenuation  in  db/cm,  c  ia  the  observed  phase  velocity, 

and  u>  is  the  angular  frequency.  These  equations  were  found  to  fit  Nolle  and 
Sleek* s  experimental  data  in  Buna-N  vulcanlzate  in  the  frequency  range  of  2  to 
10  Me  and  in  the  temperature  range  of  ~80*C  to  +80*0. 

Since  the  effective  complex  stiffness. for  longitudinal  waves  is  B*  +  4/3 u*, 

Mj  and  M2  represent,  respectively,  (B ^  +  4/3iij)  and  (B2  +  4/3u2)  when  the  longi- 
longitudlnal  velocities  and  attenuations  are  ded  in  equation  (12).  From  these 
equations  it  is  found  that  when  the  attenuation  la  low  (y  +  0)  the  wave  velocity 
equation  reduces  to  that  given  by  equations  (8)  and  (9),  viz.,  Bj  +  4/3uj  ■  ocj2. 
Therefore,  in  this  situation  the  bulk  modulus  could  be  evaluated  frem  measure¬ 
ments  on  only  the  longitudinal  and  shear  wave  velocities.  However,  for  high- 
loss  specimens,  such  as  some  of  the  nose-cone  materials  encountered  here, 
measurements  on  both  the  longitudinal  and  shear  wave  velocities  as  well  as 
attenuations  would  be  required  to  characterize  the  real  and  Imaginary  contribu¬ 
tions  to  the  bulk  modulus. 

In  addition,  measurements  on  both  the  longitudinal  and  shear  wave  velocity 
and  attenuation  Indicate  whether  longitudinal  wave  attenuation  is  due  princi¬ 
pally  to  shear  or  congressional  viscoelastic  losses.  If  the  Imaginary  components 
for  the  propagation  of  longitudinal  and  shear  waves  are  calculated  B2  +  4/3u2. 
and  u2,  respectively  ,  then  the  fraction 


■2  +  4/3w2 


(13) 
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indicates  how  much  of  the  longitudinal  wave  attenuation  is  due  to  pure  shear 
losses.  A  value  of  unity  would  Indicate  that  longitudinal  wave  attenuation  is 
due  entirely  to  shear  losses. 

Many  other  material  characteristics,  such  as  intern.*!  friction,  strain  rate 
effects,  uniformity  of  density,  existence  of  flaws,  and  specific  heats,  can  be 
obtained  by  use  of  ultrasonic  velocities  and  attenuations,  and  the  Interested 
reader  is  referred  to  the  abundant  literature  on  the  subject  (references  14,  15, 
19). 

Sound  velocities  have  many  applications  in  theoretical  and  experimental 
hydrodynamic  studies  of  impulsively  loaded  materials.  Many  of  these  applica¬ 
tions,  such  as  calculations  of  rarefaction  velocities,  momentum  transmission 
coefficients  at  impedance  mismatches,  and  pulse  lengths  of  strong  shocks,  are 
approximations  since  in  many  cases  the  corresponding  quantities  for  the  propa¬ 
gation  of  high  pressure  waves  are  not  known.  References  20  to  23  discuss,  from 
experimental  and  theoretical  points  of  view,  some  of  the  limitations  and  approx¬ 
imations  involved  in  using  ultrasonic  velocities  to  describe  the  behavior  of 
shock  waves  in  the  material  under  study.  This  report  does  not  discuss  the  con¬ 
ditions  and  limitations  of  most  of  these  approximations,  and  is  concerned  with 
some  of  the  more  quantitative  applications. 

One  method  of  studying  the  properties  of  matter  at  high  pressures  is  based 
on  the  use  of  strong  shock  waves.  A  determination  of  two  of  the  wave  parameters, 
for  example  its  velocity  of  propagation  and  the  velocity  of  the  mass  of  material 
behind  the  ahock  transition,  permits  the  pressure  and  density  in  the  shock  com¬ 
pression  to  be  found  by  conserving  momentum  and  mass  across  the  shock  front 
(Raf.  24).  The  assumptions  usually  made  are  that  a  steady  state  is  entually 
reached  behind  the  shock  front,  and  that  the  material  under  study  acts  as  a 
fluid  for  the  extremely  high  pressures  encountered  in  a  strong  shock  wave.  The 
conservation  laws  then  readily  give  a  relationship  between  the  ahock  and  parti¬ 
cle  velocities.  Also,  by  these  conservation  laws  (Rankins  -  Hugoniot  equations), 
relations  between  other  material  characteristics,  such  as  pressure  and  density, 
can  be  deduced.  The  resultant  experimental  i  arve  in  the  pressure-density  plane 
is  known  as  the  Hugoniot  curve,  which  is  defined  as  the  series  of  pressure, 
volume,  and  energy  points  generated  by  a  series  of  single  shocks  applied  to  a 
material  initially  in  a  specific  initial  state.  The  Hugoniot  equation  of  state 
for  the  material  under  study  is  then  obtained  from  several  hugoniots,  starting 
from  several  Initial  states. 
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An  Interesting  fact  that  has  emerged  from  moat  Hugoniot  measurements  is  that 
a  plot  of  the  shock  velocity,  ,  versus  particle  velocity,  up,  yields  a  straight 
line  (Ref.  24),  which  seems  to  hold  if  the  material  does  not  undergo  a  phase 
change.  The  relation  thus  experimentally  obtained 


U  -  a  +  bu_ 

■  P 


represents  the  average  behavior  of  the  large  number  of  substances  investigated. 

If  the  linear  relation  is  presumed  to  hold  all  the  way  down  to  low  pressures, 
where  the  particle  velocity  approaches  zero,  the  shock  velocity,  and  hence  the 
constant  a,  should  become  Identical  to  the  "hydrodynamic"  sound  velocity.  The 
hydrodynamic  sound  velocity,  cq.  Is  defined  as  (Ref.  22) 


2 

Co  \3p  /s 


where  PQ  is  the  pressure  corresponding  to  an  Infinitesimal  deplacement,  and  S  is 
the  entropy.  If  the  entropy  is  considereo  constant  across  the  shock  transition 
(which  is  usually  assumed),  Fawles  (Ref.  22)  indicates  that  this  is  equivalent 


/?!o\ 

\io  h 


where  H  means  that  the  derivative  is  evaluated  along  the  Hugoniot.  Since  the 
adiabatic  bulk  modulus  is  defined  as  -V  ,  equation  (16)  reduces  to 


2  m  —3 

o  p 


if  the  assumption  is  made  that  the  pressure  is  essentially  hydrostatic  in  the 
limit  of  low  pressure  and  that  the  Hugoniot  can  be  approximated  by  an  adiabat 
(see  Appendix  III).  From  the  relations  between  the  longitudinal  and  shear 
velocities  and  the  elastic  moduli  given  earlier,  equation  (17)  can  be  expressed 
as  (see  Appendix  III) 

co  "  C1  '  4,3  ct 


(18) 
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Therefore,  the  extrapolation  point  a  ■  cq  can  1 2  estimated  from  hydrostatic  P-V 
data  corrected  to  constant  entropy  (Ref,  24),  or  from  the  longitudinal  and  trans¬ 
verse  velocities  in  the  materials.  The  comparison  between  calculated  values 
using  equation  (18)  and  extrapolated  values  is  good  for  most  metals,  and  some 
values  are  given  in  Appendix  III,  As  shown  in  equation  (18),  the  calculated 
value,  cq,  is  lower  than  the  longitudinal  velocity  in  the  material,  and  this  is 
generally  true  of  the  extrapolated  values,  if  the  Hugoniot  experiments  are  per¬ 
formed  at  significantly  high  pressures.  Figure  34  Illustrates  the  typical  vari¬ 
ation  between  the  longitudinal  velocity  and  the  extrapolated  value  of  shock 
velocity. 

If  the  linear  relation  (equation  14)  is  considered  as  an  expansion  of  the 
shock  velocity,  Paul  and  Varschauer  (24)  show  that  the  slope,  b,  should  be  given 

as 


which  can  be  further  reduced  to 

b-.M  [l  +  (£)] 


(19) 


(20) 


The  quantity  ^  -JL J  ,  which  is  the  change  of  the  adiabatic  bulk  modulus  with 
pressure,  can  be  obtained  from  the  pressure  dependence  of  the  longitudinal  and 
transverse  velocities.  Table  V  shows  the  agreement  between  the  values  of  the 
slope  calculated  fro a  ultrasonic  data  compared  to  data  obtained  from  shock  ex¬ 
periments  as  listed  by  Paul  (Isf.  24). 

The  data  obtained  from  shock  measurements  is  generally  expressed  In  equation 
f  state  fora  as  an  expansion  in  tits  compression  (Safe.  23,26) 


(21) 
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Table  V 

COMPARISON  OF  THE  SLOPES  OF  THE 
SHOCK  VELOCITY  VERSUS  PARTICLE  VELOCITY  CURVE 


Material 

Ultrasonic,  b 

Shock,  b 

NaCl 

1.35  +  .05 

1.43  +  .02 

Aluminum 

1.54  ±  .05 

1.2  to  1.6 

Copper 

1.64 

1.50 

Silver 

1.79 

1.59 

Gold 

1,85 

1.56 

Sodium 

1.15 

1.2 

* 

$ 
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P 

where  y*  -  ~  -1  is  the  compression,  and  the  zero  subscript  indicates  that  the 

derivatives  are  to  be  evaluated  at  the  initial  conditions.  Since  many  experi¬ 
ments  are  made  with  pressures  on  the  order  of  tens  of  kilobars  or  greater,  PQ  is 
usually  neglected.  Again,  the  first  derivative  can  be  related  to  the  bulk  modu¬ 
lus,  which  gives 


(22) 


which  can  be  further  related  to  the  adiabatic  velocities 


-  PCc l 


(23) 


(For  the  assumptions  involved  in  this  derivation,  see  Appendix  III).  The  longi¬ 
tudinal  and  transverse  velocities  can  then  be  used  to  obtain  a  theoretical 
constant  useful  in  describing  experimental  hydrodynamic  data.  The  temperature 
dependence  of  the  longitudinal  and  transverse  velocities  allows  a  convenient 
calculation  of  a  constant  to  be  used  in  experiments  under  varying  temperature 
environments. 

Normally,  the  equation  of  state  used  to  describe  impulsively  loaded  materials 
contains  an  additional  term  to  equation  (21)  which  involves  the  pressure  depend¬ 
ence  of  the  thermal  energy,  E,  of  the  solid 


P  •  ay*  +  by’2  +  cy’3  + 

po 


(24) 


where  y  is  obtained  from  the  Mie-Grueneisen  equation  of  state  (Ref.  22).  This 
parameter  enters  into  all  calculations  concerning  the  description  of  shock  waves 
propagating  through  materials.  As  shown  by  Paul  (Ref.  24),  the  Grueneisen  ratio 
is  defined  as 


Y 


(25) 
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1 


j 

I 

Using  thermodynamic  relations  and  the  relations  among  the  velocities f  this  ratio  j 
can  be  expressed  for  isotropic  materials  as  (see  Appendix  IV)  l 


i 


where  8  is  the  volume  coefficient  of  expansion,  and  cp  is  the  specific  heat.  The 
coefficient  calculated  in  this  way  for  metals  compares  closely  with  that  calcu¬ 
lated  from  experimental  shock  wave  data  (Table  VII,  Appendix  IV). 

In  another  application,  composite  materials  are  used  to  reduce  the  magnitude 
of  a  high  pressure  Impulse  at  a  solid  interface.  This  necessitates  a  knowledge 
of  the  reflection  and  transmission  coefficients  of  a  stress  wave  at  the  boundary. 
Allen  (Ref,  27)  discusses  the  applicability  of  a  nonlinear  approach  to  the  calcu¬ 
lation  of  the  transmission  coefficient,  which  in  general  exceeds  that  calculated 
from  the  linearised  equations  of  motion.  However,  the  derivation  by  the  linear¬ 
ized  approximation  gives  the  transmission  (t)  and  reflection  (r)  coefficients  in 
terms  of  the  acoustic  Impedances  of  the  two  solids  as  (Ref.  24) 


r  _  Po2C2  ~  P01C1 

h  - 

*1 

"ol'l  +  ”o2C2 

z1  + 

Z2 

(27) 

2p0l=I 

2Zt 

'  ‘  polcl  +  "ozS 

-  — 

h 

where  the  wave  is  propagating  from  medium  1  into  medium  2,  and  the  c1  and  c2  are 
the  shock  velocities  in  the  two  media.  Longitudinal  velocities  offer  convenient 
approximations  in  the  calculations  of  transmission  coefficients  when  the  equa¬ 
tions  of  state  are  not  known  and  can  be  used  as  approximations  in  equations  (27). 

The  data  presented  earlier  on  the  longitudinal  velocity  through  laminated 
materials  as  a  function  of  the  angle  of  Incidence  also  offers  an  attractive 
approach  to  the  problem  of  shock  reduction,  since,  as  previously  shown,  the 
velocity  is  a  sensitive  function  of  the  incident  angle  for  some  layered  materials. 
A  shock  wave  could  then  be  attenuated  by  a  material  with  a  continuously  varying 
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lamination  angle,  thus  eliminating  excess  bonds  between  materials  and  possibly 
preserving  desirable  ablation  characteristics.  However,  the  approach  ie  still 
tentative,  and  some  theoretical  work  is  needed  to  correlate  the  present  acoustic 
measurements  to  the  actual  condition  of  high-intensity  shock  propagation. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  temperature  dependence  of  longitudinal  and  transverse  velocities  offers 
a  convenient  means  of  calculating  many  material  properties  as  a  function  of 
temperature.  The  knowledge  of  mechanical  properties  at  high  temperatures  is 
becoming  an  Increasingly  important  requirement  in  hydrodynamic  studies,  and  such 
knowledge  is  enhanced  as  experimental  techniques  and  theoretical  approaches  be¬ 
come  more  sophisticated. 

Many  of  the  results  obtained  from  sound  velocities  can  be  duplicated  by 
other  means,  such  as  static  tests  in  the  case  of  the  elastic  moduli.  However, 
ultrasonic  measurements  offer  a  check  on  the  static  results  and  perhaps  apply 
more  directly  to  hydrodynamic  experiments  because  of  the  dynamic  character  of 
both. 

Another  consideration  is  that,  in  plastics  relaxation  effects  can  be  lm> 
portant  so  that  no  simple  relationship  exists  between  the  various  median-  al 
properties  of  the  material  as  determined  by  static  and  dynamic  measurements. 

That  is,  mechanical  properties  such  as  the  bulk  modulus  are  a  function  of  the 
strain  rate.  For  example,  the  adiabatic  bulk  modulus  calculated  from  velocity 
measurements  can  be  considered  to  consist  of  two  parts!  one  associated  with  an 
infinite  relaxation  time,  which  is  a  function  of  the  lattice  spacing,  and  one 
given  by  a  finite  relaxation  time,  relatin,  to  structural  rearrangements.  Such 
effects  may  be  Important  in  understanding  the  relationship  between  sound  veloci¬ 
ty  measurements  as  obtained  from  the  Hugoniot  equation  of  state  and  those  from 
ultrasonic  measurements.  Sound  velocities  at  different  frequencies  should  be 
useful  in  studying  these  relaxation  phenomena. 

Other  parameters,  such  as  the  Grusnelsen  ratios  calculated  from  sound 
measurements  and  shock  measurements  can  be  compared  to  give  further  theoretical 
insight  of  the  hydrodynamic  modal  of  the  material. 

As  mentioned  earlier  (aquation  20) ,  the  slope  of  the  shock  velocity  - 
particle  velocity  curve  for  metals  can  be  estimated  quite  closely  by  means  of 
ultrasonic  pressure  measurements .  Similar  ultrasonic  experiments  should  be 
performed  for  the  increasingly  important  plastics  and  nose- cone  materials  studied 
in  high  i^>ulse  experiments.  Also,  with  the  proper  experiment  and  correct  inter¬ 
pretation  of  data,  the  pressure  dependence  of  ultrasonic  velocities  might  be 
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used  to  explain  the  discrepancy  In  the  lover  pressure  region  of  the  shock  veloc¬ 
ity/particle  velocity  curve  as  shovn  in  figure  34. 

Theoretical  investigations  should  be  organized  in  two  particular  directions. 
First,  the  relations  between  the  adiabatic  longitudinal  velocities  and  the 
mechanical  properties  of  the  laminated  materials  should  be  Investigated  in  de¬ 
tail  to  correlate  experimental  data  to  hydrodynamic  parameters  such  as  the 
Gruenelsen  ratio,  reflection  coefficients,  etc.  Secondly,  the  temperature  de¬ 
pendence  of  longitudinal  velocities  should  be  thoroughly  investigated  as  to 
application  in  theoretical  predictions  of  high-pressure  experiments.  For  exam¬ 
ple,  the  results  of  the  present  experiments  show  that,  in  general,  the  longi¬ 
tudinal  velocities  in  most  nose-cone  materials  decrease  with  increasing  tempera¬ 
ture.  This  Indicates  that  rarefaction  velocities  would  likewise  decrease  with 

•» 

increasing  temperature  (the  internal  energy  and  temperature  are  higher  behind 
the  shock  front)  and  consequently  affect  such  parameters  as  shock  width,  since 
the  rarefaction  wave  effectively  terminates  the  impulse  loading.  (The  impulse 
is  determined  by  the  relative  wave  front  and  rarefaction  velocities.  For  a 
temperature  gradient  this  relation  changes  if  a  velocity  gradient  exists  for 
the  rarefaction  wave.) 

Future  ultras  ooic/rraterials  investigations  at  AFWL  will  involve  a  measure¬ 
ment  of  shear  velocities  as  a  function  of  both  temperature  and  pressure  for  the 
materials  reported  here.  Longitudinal  measurements  will  likewise  be  made  as  a 
function  of  uniaxial  stress  and  will  be  extended  in  temperature  in  an  attempt 
to  shed  light  an  some  of  the  above  problem  areas. 
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APPENDIX  I 

INTERNAL  TEMPERATURES  OF  THE  SAMPLES 

During  the  course  of  the  Investigation,  sample  temperatures  were  controlled 
by  either  a  thermostatically  controlled  oven  or  an  oil  bath.  The  teiq>erature 
was  determined  by  reading  the  temperature  of  the  oil  or  the  air  temperature  near 
the  sample  and  not  the  Internal  temperature  of  the  sample  Itself.  Therefore,  to 
eliminate  any  temperature  gradients  in  the  sample  it  was  necessary  to  determine 
the  time  required  for  the  sample  to  obtain  a  uniform  temperature.  An  experi¬ 
mental  study  was  conducted  to  determine  the  internal  temperatures  of  the  samples 
as  a  function  of  time. 

Samples  were  prepared  by  drilling  small  holes  in  various  places  of  the 
sample  geometry  and  by  inserting  iron- cons tantan  thermocuple  wires  as  shown  in 
figure  35.  Internal  temperature  gradients  were  then  determined  by  measuring  the 
output  of  each  thermocouple.  The  measurements  were  made  by  first  bringing  the 
environment  to  the  desired  temperature.  Inserting  the  sample,  and  measuring  the 
temperature  at  each  position  in  the  sample  as  a  function  of  time. 

A  variety  of  temperatures  were  used  as  the  ambient  temperature,  as  indicated 
in  the  following  graphs,  although  it  was  observed  that  for  the  range  studied  the 
time  required  for  equilibration  was  fairly  Independent  of  initial  temperature. 

Figures  36  and  37  show  the  temperature  curves  obtained  by  heating  samples  c? 
Plexiglass  and  low-density  polyethylene  in  a  70*C  oil  bath.  It  is  observed  that 
temperature  equilibrium  is  readied  in  about  one  hour  In  both  materials. 

Figures  38  and  39  show  die  internal  temperature  variation  in  AVCO  fiberglass 
and  OTWR  at  temperatures  of  217*C  and  2i2*C,  respectively.  For  these  measure¬ 
ments  die  oven,  rather  than  the  oil  bath,  was  used  to  control  the  temperature 
environment.  It  le  noted  that  the  same  form  of  temperature  curve  exists  for 
both  methods,  which  implies  that  the  waiting  period  is  Independent  of  the  origi¬ 
nal  temperature  over  this  range  of  temperatures. 


77 


AfVL- Tit-65-168 


3.  3UniVU3dMli 


Flgur*  36.  Internal  Taaparatxiraa  of  Plexiglass  Saaple;  Initial  Taaparatura  of  Bach. 
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r Inure  37.  Internal  Ta«pe rat urea  of  Polya thy lana  Samplaa;  Initial  Tanparatura  of  Bath, 
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Internal  Teape returns  of  AVCO  Fiber  Glass;  Initial  Tenperature  of  Oven 

217*  ♦  *C 


•2 


212*  +  1  *C 
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APPENDIX  II 

DERIVATION  OF  EQUATIONS  OF  MOTION  FOR  ULTRASONIC  WAVES 

The  generalized  form  of  Hook's  law  up  to  the  elastic  limit  is  given  by 
(Ref.  l4). 

\ 

°xx  "  -ll£xx  +  c12£yy  +  c13ezz  +  ci4Eyz  +  cisezx  +  c16exy 

0  -  c,.e  +  c,,e  +  c„„e  +  cn  e  +  c._e  +  c  e 

yy  21  xx  22  yy  23  zz  24  yz  25  zx  26  xy 

a  _  c  £  +  c,  e  +  c.  e  +  c  ,e  +  e„„e  +  c,  e  (281 

zz  31  XX  32  yy  33  zz  3 4  yz  35  zx  36  xy  K  ' 


0  ■  c.,e  +  c,  e  +  c  e  +_c-r  ,e  +  c,  e  +  e. 

yz  41  xx  42  yy  43  zz"^h4  yz  45  zx  46  xy 


a  *  c,  e  +c.e  +  c_  £  +ce  +ce  +  cCFe 
zx  51  xx  5?  yy  53  zz  54  yz  55  zx  56  xy 


o  -  c  e  +  c,  e  +  c,,e  +  c„  e  +  c  e  +  c  £ 
xy  61  x..  62  yy  63  zz  64  yz  65  zx  66  xy 


where  the  coefficients  are  the  elastic  constants  of  the  material.  The  first 


letter  in  the  subscript  of  the  stress  ^0ijj  denotes  the  direction  of  the  stress 
and  the  second  letter  defines  the  plane  in  which  it  is  acting.  The  identical 
subscripts  on  the  strain  ^eij^  refer  to  extensional  and  shearing  strains.  For 
the  elastic  energy  to  be  an  mivalued  function  of  the  strain*  any  coefficient 
c^j  must  equal  the  coefficient  Cjj  (Ref.  14).  The  number  of  independent  co¬ 
efficients  is  thus  reduced  from  j6  to  21.  If  the  material  has  axes  or  planes 
of  symmetry,  the  required  number  of  independent  coefficients  is  correspondingly 
reduced.  For  •'xample,  cubic  symmetry  requires  only  three  independent  constants* 

c  .  c  ,  and  c  . 

11*  ] 2 1  14 

In  an  isotropic  material  the  values  of  the  coefficients  must  be  independent 
of  the  chosen  coordinate  system.  This  constrains  the  matrix  to  two  coefficients 
which  ate  commonly  known  as  the  Lame  constants,  X  and  y.  Kolsky  Qtef  14)  than 
shows  that 
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C  m  C  M  C  *  r  ■  r  s  1 

12  13  C2 1  C2  3  C31 


c, ,«  c,  *  c  ■  u 
**•*  55  66 


C11  "  c22  *  c3  3  *  ^  +  2u 


Equations  (28)  may  then  be  written  as 


a _ -  XA  +  2ye _ ;  o  -  XA  +  2ye  ;  a  ■  XA  ■+  2ye 

xx  xx  yy  yy  zz  zz 


O  '  «  ye  ;  or  ■  ye  ;  o  «  ye 

ya  yz  zx  zx  xy  xy 


where  A  «  e _  +  e  +  e  represents  the  change  in  volume  and  is  called  the 

xx  yy  zz 

dilatation. 

To  obtain  the  equation  of  motion  for  an  elastic  medium  the  net  forces  acting 
on  a  small  parallelepiped  of  volume  6x<Sy6z  are  equated,  to  the  rate  of  change 
of  the  momentum  in  any  direction.  Consider  the  net  force  in  the  x  direction. 
This  is  given  from  the  definition  of  the  subscripts  on  the  stresses  as 


Vx.  *  -  axx6*Sz  +  (v  +  if4*)' 

I  9°™  \ 

-  o  6xAz  +  lo  +  -r— 6zJ6x6y  -  or  6x6y 
xy  \  xic  oz  /  xz 


|6x6y  -  or  4*«y 


where  the  origin  is  placed  at  one  corner  of  the  parallelpiped.  By  Newton's 
second  law  of  motion  the  net  force  is  equated  to  the  rate  of  change  of  mass 
times  velocity  in  the  x-dir'  ction 


(if  +  if  +  if)4,,4>,4x  • 
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where  p  is  the  density  and  u  is  the  displacement  in  the  x  -  direction.  Similar 
equations  apply  when  the  displacements  in  the  y  and  z  -  directions  (v  and  v, 
respectively)  are  considered.  The  resulting  equations  for  the  three  directions 
then  become 


P 


3t2 


3  o 

—22 

3y 


,2  3a  ,  3a  .  3 a 
m  __22  +  _JOL  +  .  .  ja 


3t2 


3x 


3y 


3z 


ii? 

3t2 


3o  3o  30 
3x  3x  3a 


(33) 


Substituting  from  equation  (30)  gives  for  the  equation  of  motion  in  the  x  - 
direction. 


(34) 


Now  consider  a  point  in  the  undisplaced  position  which  has  coordinates  (x  +  fix, 
y  +  6y,  z  +  fiz),  and  let  the  displacement  which  it  lias  undergone  have  components 
(u  +  6u,  v  +  Sy,  w  +  6w) .  If  fix,  5y,  and  5z  are  sufficiently  small, 

.  3u-  .  3u„  .  3u„_. 

6u  .  ~6x  +  -^fiy  +  -^fiz 


i  3v.  ,  3v-  .  3vf_ 

fiv  -  —fix  +  -g^y  + 


6W  "  f 6X  +  l?6y  +  It*1 


where 


3u  3v  3w 


represent  the  fractional  expansions  and  contractions  of 


3x*  3y*  3s 

infinitesimal  line  elements  parallel  to  the  x,  y,  z  direction,  respectively. 


and  are  respectively  denoted  by  e 


•  € 

xx  yy' 


The  JSL  x  jiu  .  3w, 

„„  'zz*  P  8  3y  3z*  3-;  +  3x 

•—£  +  correspond  to  the  components  of  shear  strain  in  the  planes  denoted  by 
their  suffixes.  These  quantities  are  defined  as 
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e  -  +  i2 

yz  3y  3z 


3u  +  3w 
'zx  3z  3y 


(36) 


_  3v  .  3u 

xy  3x  3y 


Rotations  of  the  element  about  a  point  are  given  as  the  difference  of  the 
partial  derivatives,  whi^h  Kolsky  (Ref.  14)  defines  as 


2w  *» 
z 


3w 

3v 

3y  " 

3z 

3v 

3w 

3z  ~ 

3x 

3v 

3u 

3x  “ 

3y 

(37) 


Using  these  definitions  the  equations  of  motion  become 

-  (X  +  y)f  +  y?2u 

3t2  3x 


9~-  -  (X  +  y)~  +  yV2v 
3t2  * 


(38) 


pl^H  ■  (X  +  y)|~>  +  yy^ 

3t2  38 
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Manipulation  of  these  equations  gives 


—  -  (X  +  2y)72A 
3t2 


(39) 


The  first  equation  shows  that  the  dilation  A  is  propagated  through  the  medium 

1/2 

,  while  the  second  corresponds  to  a  rotation  propa- 
1/2 

waves  in  the  interior  of  an  Isotropic 


with  velocity 


h- 


hi 


gated  with  a  velocity 


Mi/z 

.  Thus, 


elastic  solid  may  be  propagated  with  two  different  velocities.  Waves  involving 

1/2 


no  rotation  travel  with  velocity  c^ 
dilation  travel  with  velocity  cfc  - 


and  waves  involving  no 


1/2 
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APPENDIX  III 

COMPARISON  OF  ADIABATIC  LONGITUDINAL  VELOCITY 
AND  THE  EXTRAPOLATED  SHOCK  VELOCITY 


As  shown  in  Appendix  II  the  longitudinal  wave  velocity  for  small  stresses 
in  a  homogeneous ,  unbounded,  isotropic  medium  is  given  by 


(40) 


where  A  and  u  are  the  adiabatic  Lamife  constants.  As  previously  shown,  for 
infinitesimal  displacements  the  equations  governing  a  longitudinal  wave  propa¬ 


gating  in  the  x  -  direction  (e 


yy 


zz 


0)  are  given  as 


o  -  (A  +  2u)e 

XX  XX 


(41) 


Since  e  »  -~t  where  u  is  the  displacement  of  a  point  in  the  body  by  the  stress 

XX  dX 

o^,  Morgan  (Ref.  26)  shows  that  equation  (40)  can  be  written  as 


l  VP  3u  *  9p 


(42) 


This  equation  applies  only  when  the  induced  stresses  imposed  on  the  sample  by 
the  wave  are  small. 

If  a  large  stress  is  induced  in  the  material  by  some  means  such  as  that 


produced  by  a  high  velocity  impact,  the  wave  velocity,  UR  ■M- wiu  b* 
constant  and  independent  of  the  stress  level  if  the  induced  stress  is  less  than 
the  elastic  limit.  In  this  region  the  elastic  constants  alone  determine  the 
slope  of  the  one-dimensional  strain  curve  which  has  a  slope  of  A  4-  2y.  When 
plastic  yielding  occurs  the  curve  tends  to  parallel  the  hydrostat  which  has  a 
slope  of  ^A  (Ref.  28) .  This  transition  regio.  where  the  strain  slope 

changes  from  ^A  +  is  usually  called  the  Hugoniot  elastic  limit. 
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In  this  discussion  it  is  assumed  that  shear  effects  have  an  important  effect 
on  the  equation  of  state.  Consider  the  stress-strain  plot  in  figure  AO.  Below 
the  elastic  limit  (point  A)  any  stress  disturbance  will  propagate  with  a  constant 
velocity  as  normally  measured  in  the  laboratory.  For  pressures  between  A  and  B, 
where  B  refers  to  the  inflection  point  at  which  ^  starts  increasing  with  in¬ 
creasing  p,  the  stress  wave  propagation  velocity  decreases  with  increasing  stress. 
The  wave  will  thus  break  up  into  a  two  wave  structure  -  the  elastic  portion 
traveling  at  the  dilatational  velocity  and  the  plastic  portion  traveling  with  a 


-GO 
’•&) 


at  that  pressure,  where  refers  to  the 


velocity  given  by  the  value 
one  dimensional  stress. 

From  the  point  B  to  C,  |-^|  increases  with  increasing  stress  level,  with 
the  net  effect  of  an  unstable  shock  wave  preceeded  by  an  elastic  and  plastic 
precursor.  For  pressures  greater  than  that  corresponding  to  point  C  the  veloci¬ 
ty  of  propagation  of  the  disturbance  is  greater  than  the  elastic  wave  velocity 
and  a  single  shock  front  occurs. 

For  a  fluid  or  material  in  which  p  is  xero  the  stress  strain  curve  of  figure 
AO  does  not  show  the  Inflection  behavior  and  is  a  monotonlcally  increasing 
function  of  the  strain  after  the  elastic  limit  is  exceeded  (if  no  phase  transi¬ 
tions  are  present).  Morgan  (Ref.  26)  states  that  in  this  region  the  bulk  modu¬ 
lus  it.  then  given  by  1  rather  than  x  *&■)• 

Morgan  (Ref.  26)  discusses  the  extrapolated  shock  velocity  from  a  theoreti¬ 
cal  point  of  view  for  two  situations.  Bis  predictions  indicate  that  for  measure¬ 
ments  taken  in  the  region  where  the  inflection  due  to  shear  rigidity  occurs 


(see  figure  AO) ,  the  use  of  m.  as  the  wave  velocity  and  the  extrapola¬ 
tion  of  the  data  to  xero  pressure' yields  the  sound  velocity  as  normally  measured 
for  very  small  stresses.  The  extrepolated  velocity  should  then  be  given  as 


».<»  -  „>  -  c,  -VFH, .  o  “3) 

The  aVCO  report  (Ref.  26)  states  that  for  pressures  well  above  the  inflec¬ 
tion  point  the  shear  rigidity,  p,  can  be  neglected  relative  to  the  magnitude 
the  volume  compressibility,  X,  as  previously  noted.  In  this  region  it  has  been 
found  experimentally  that  the  shock  velocity  is  usually  a  linear  function  of 
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Figure  40.  Dynamic  Stress-Strain  Curve  for  a  Material  with  Shear  Rigidity 
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the  particle  velocity.  Thus*  using  this  hypothesis,  the  report  shows  that  the 
extrapolated  value  of  shock  velocity  should  be 


Ug(P 


(44) 


From  the  relations  given  in  equation  (39),  equation  (44)  can  be  expressed  as  a 
function  of  longitudinal  and  transverse  velocities  at  atmospheric  pressure. 
Table  VI  lists  some  experimentally  extrapolated  velocities  Ug(P  *  0)  and  the 
corresponding  predicted  values,  c*,  of  equation  (44)  as  listed  in  die  AVCO 
report.  It  is  noted  that  AVCO 'a  predicted  values  are  generally  lower  than  the 
actual  extrapolated  values,  particularly  in  the  case  of  Beryllium. 

The  table  also  lists  the  adiabatic  longitudinal  velocities  at  atmospheric 
pressure,  c#,  and  some  values  under  the  heading  c  which,  will  presently  be 
explained. 


Table  VI 

COMPARISON  OF  EXTRAPOLATED  SHOCK  VELOCITIES  AND  CALCULATED  VALUES 


ci 

U#(P  -  0) 

c* 

co 

Metal 

a/jis*£ 

mn/usec 

mm/usec 

ssLasi 

Gold 

3.24 

3.059 

2.77 

2.94 

Beryllium 

12.89 

7.975 

2.96 

7.80 

Hagnsalue 

5.77 

4.493 

3.82 

4.57 

nickel 

6.04 

4.667 

4.32 

4.95 

Tin 

3.32 

2.668 

2.39 

2.64 

Zinc 

4.21 

3.042 

2.43 

3.15 

Aluminum 

6.42 

5.32 

4.75 

5.37 

Silver 

3.65 

3.215 

2.85 

3.13 
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The  predicted  value  of  the  extrapolated  shock  velocity  as  given  by  equation 
(18)  can  be  derived  using  the  assumption  of  the  Hugoniot  equation  of  state  that 
the  pressure  can  be  expressed 


(45) 


> 


neglecting  the  term  involving  the  internal  energy.  Ass«;alng  that  P(y')  is 

continuous,  that  is,  that  inflective  or  transitive  behavior  does  not  exist  the 

definition  of  the  shock  velocity  can  be  used  to  calculate  the  value  in  the 

limit  of  ?•*?,»'♦  o 
o 


(46) 


Using  the  relation 


equation  (46)  becomes 


(47) 


(48)  I 

i 
s 


where  the  derivative  is  to  be  evaluated  at  the  initial  conditions.  From  the 
definition  of  the  adiabatic  bulk  modulus,  equation  (48)  can  be  expressed  as 

B 

U2(P  -  0)  -  —  (49) 

s  po 
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This  last  step  results  In  the  sane  extrapolated  value  of  shock  velocity  that  Is 
predicted  In  the  AVCO  report  (equation  44)  if  It  is  assumed  that  the  bulk  modu¬ 
lus  at  high  pressure  la  given  by  X  and  should  consequently  extrapolate  to 
Xp  m  under  normal  conditions.  It  is  assumed  here  that  the  extrapolated  value 
corresponds  to  a  hydrostatic  compression  and  thus  is  related  to  the  bulk  modulus 


X  +  "vj  as  measured  under  atmospheric  pressure.  This  assumption  is  partially 
supported  by  Bridgman's  measurements  on  the  isothermal  compressibilities.  His 
data  in  the  upper  pressure  range  generally  agree  to  within  a  few  percent  of  the 
data  in  the  lower  pressure  range  (^100  kb)  of  dynamic  experiments  (for  example, 
see  Walsh's  comparison  (Ref.  29)  and  some  of  Bridgman's  original  data  (Ref.  30)) . 


Bridgman's  data  for  the  compression  ~  are  fitted  by  a  quadratic  expression 

o 

A  V  q 

in  pressure:  -  -aP-bPz,  where  a  is  the  isothermal  compressibility  (reciprocal 

vo  f  2  1  . 

of  isothermal  bulk  modulus  |X  +  —  y  and  b  is  a  constant.  If  Bridgman  s  data, 
corrected  to  the  adiabat  (to  correspond  to  the  approximately  adiabatic  dynamic 
measurements),  are  assumed  to  be  9  continuous  extension  of  the  dynamic  measure¬ 
ments,  than  the  slope  of  static  curve  near  zero  pressure  should  correspond  to 
the  extrapolated  slope  of  the  dynamic  curve.  Under  these  assumptions  the  extra¬ 
polated  shock  velocity  should  be  given  as  the  square  root  of  (X  +  2/3y)p  _  e, 
which  by  equation  (9)  can  be  further  related  to  longitudinal  and  shear  veloci¬ 
ties. 


U2(P  -  o)  «  c2  -  4/3  c 2  (50) 

S  l  * 

which  was  previously  defined  as  cQ.  The  assumptions  used  in  evaluating  this 
derivative  near  zero  pressure  also  apply  to  the  justification  of  equation  (23) 
since  the  derivative  in  both  cases  is  the  same. 

The  actual  situation  is  much  more  complicated  than  this  simple  treatment, 
and  some  of  the  assumptions  used  here  are  certainly  not  valid  in  all  cases. 
However,  it  is  interesting  to  note  in  Table  VI  the  fair  agreement  between  the 
predicted  values  of  extrapolated  shock  velocity,  cQi  from  equation  (50)  and 
those  obtained  by  experimental  extrapol  ation  U8(P  -  o) . 

Equation  (50  rests  upon  the  assumptions  that  (1)  the  data  can  be  expressed 
by  s  smooth  curve  of  the  form  of  equation  (45),  (2)  that  the  Hugonlot  is 
sufficiently  near  tha  adiabat,  so  that  U2  «  ^  *****  00 

clone  or  inflections  occur  in  the  observed  pressure  range,  (4)  that  the  shock 


93 


AJWL-TR-65-188 


velocity  is  given  by  a  linear  relation  with  the  particle  velocity,  and  (5)  that 
the  internal  energy  tern  in  the  pressure  can  be  neglected. 

Assumptions  (1)  and  (3)  are  in  general  justified  for  pressure  experiments 
above  100  kb  (Ref.  29) .  The  Hugoniot  is  usually  within  a  few  percent  of  the 
adiabat  so  that  assumption  (2)  is  a  fair  approximation.  As  mentioned  earlier, 
the  linear  relationship  between  shock  and  particle  velocity  represents  only  the 
average  behavior  of  a  large  number  of  substances.  Therefore,  .umptlon  (4) 
must  be  evaluated  for  the  individual  experiment.  The  effect  at  the  energy 
term  has  on  the  predictions  of  equation  (SO)  is  not  new  known,  but  might  be 
responsible  for  some  of  the  discrepancies  in  Table  VI. 
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APPENDIX  IV 

EVALUATION  OF  GRUENEISEN  RATIO  FROM  A  KNOWLEDCE 
OF  SOUND  VELOCITIES 


The  Grueueisen  constant  is  defined  as  (Ref.  26), 


where  V  is  the  specific  volume  —  ,  E  is  the  internal  energy,  and  is  the 
specific  heat  at  constant  volume.  Now, 


(52) 


where  8  is  the  volume  coefficient  of  expulsion  and  ia  the  Isothermal  bulk 
modulus.  It  is  desired  to  express  equation  (52)  _’n  terms  of  the  acoustic  veloc¬ 
ities,  which  are  obtained  under  adiabatic  conditions.  This  implies  that  the 
associated  elastic  moduli  must  be  the  adiabatic  rather  than  the  IsoCieimal 
moduli.  With  this  in  mind  and  applying  the  relation  between  the  adiabatic  and 
isothermal  bulk  modulus  the  last  part  of  equation  (52)  becomes. 


B 

s 


c 

v 


(53) 


Using  this  relation  and  the  first  part  of  equation  (52)  yields  for  equation  (51) 


y(V,T) 


PCT 


(54) 
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The  adiabatic  bulk  modulus  is  defined  in  terms  of  the  Lairife  constants  for  an 
isotropic  homogeneous  medium  (Ref.  1)  as 


B 

s 


(55) 


The  dilatational,  c£,  and  transverse,  velocities  in  an  isotropic  medium  are 
shown  in  Appendix  II  to  be 


c 


% 


J  ' 

so  that  with  equations  (55)  and  (56)  equation  (54)  becomes, 


y(V,T)*-|-  -  4/3  c* 


(56) 


(57) 


Equation  (57)  gives  the  value  of  the  Grueneisen  constant  for  isotropic  materials 
as  determined  from  sound  velocity  measutements  and  compares  closely  with  that 
computed  by  other  means.  The  following  table  shows  some  values  of  the  ratio  y^ 
computed  from  shock  measurements  and  extracted  from  the  AVCO  report  (Ref.  26)  and 
the  corresponding  values  yy  calculated  from  sonic  measurements.  The  necessary 
quantities  used  in  equation  (57)  were  handbook  values. 
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TABLE  VII 

COMPARISON  OF  GRUENEISEN  RATIOS  CALCULATED  FROM  SHOCK  MEASUREMENTS 
TO  THOSE  CALCULATED  FROM  ULTRASONIC  VELOCITIES 


Material 

Is 

In 

Beryllium 

1.17 

1.29 

Copper 

2,04 

2.40 

Magnesium 

1.46 

1.56 

Tin 

2.03 

2.05 

Zinc 

2.38 

2.30 

24  ST  Aluminum 

2.13 

2.25 

Sliver 

2.47 

2.45 
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APPENDIX  V 

CALCULATION  OF  TRANSMISSION  COEFFICIENTS 
AT  IMPEDANCE  MISMATCHES 


The  purpose  of  this  section  will  be  to  derive  the  correction  factor  C  used 
in  equation  (4) .  Consider  the  following  sketch: 


SAMPLE 


BUFFER 


_js  A*!  _ 

aT 

Zb 

A|2 

Z, 

Af“ 

Zr 

7 

V 

BUFFER 


Here  a  sample  of  Impedance  Z  (»  p  c  )  is  inserted  between  two  buffers  of  im- 

pedance  Z^(»pBCg) .  Neglecting  the  small  loss  through  the  buffers  the  input 

amplitude  A  and  the  output  A  is  measured.  From  the  ratio  of  these  two  signals 
o  t2 

it  is  desired  to  calculate  the  amplitudes  A  and  A  '  which  give  a  measure  of 

C1  C1 

attenuation  of  the  signal  through  the  specimen.  2g 

At  surface  S,  the  transmission  ratio  is  given  by  t  -  - — so  that  the 

h  +  ZS 

amplitude  A  is  given  by 
C1 

2ZS 

At  "  rrr  Ao  (58> 

1  B  S 


98 


AFWL-TB.-65-188 


and  at  surfaces  S2, 


A  is  given  by 
z2 


(59) 


Now  the  natural  logarithm  of  the  two  amplitudes  A'  /A  is  by  definition,  the 

tj 

attenuation  through  the  sample  in  Nepers,  -  a  where  e°  is  the  ratio  of  the  am¬ 
plitude  at  a  point  x  to  that  at  the  point  x  +  dx.  The  corresponding  attenuation 
in  decibels  per  unit  length  L  is  L  -  201og10(ea)  «■  8.686  o.  By  definition 


where  d  is  the  sample  thickness. 
Taking  logs 


ad  log1Qe  -  log 


10  A 


1 

20alogl0e  -  -  201og10  jr*- 


(60) 


(61) 


or 


L 


~  201og 


1 


Substituting  in  the  values  of  A 

cl 


and  A  '  in  terms  of  known  quantities  A  , 


gives  for  L 


(62) 
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This  correction  factor  C  is  the  one  defined  in  equation  (4)  and  accounts  for 
the  impedance  mismatch  at  the  boundaries .  For  materials  in  which  Idle  impedances 
are  similar  to  that  of  aluminum  the  correction  factor  is  negligible.  However, 
for  the  foam  materials  the  impedance  mismatch  between  the  buffer  and  specimen 
was  very  large  so  that  the  correction  factor  was  on  the  order  of  20  db. 

This  correction  does  not  account  for  other  influences  such  as  surface  de¬ 
fects,  impedance  f  the  bond,  etc.  which  must  be  evaluated  in  other  ways.  How¬ 
ever,  for  the  high  loss  epoxy  foam  discussed  in  Section  III,  the  term  C  in 
equation  (63)  was  found  to  be  the  most  significant  factor  in  correcting  for  the 
total  loss. 
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